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ABSTRACT 

The damping of the harmonic oscillator is studied in the framework of the Lindblad 
theory for open quantum systems. A generalization of the fundamental constraints on 
quantum mechanical diffusion coefficients which appear in the master equation for the 
damped quantum oscillator is presented; the Schrodinger, Heisenberg and Weyl-Wigner- 
Moyal representations of the Lindblad equation are given explicitly. On the basis of these 
representations it is shown that various master equations for the damped quantum oscilla- 
tor used in the literature are particular cases of the Lindblad equation and that not all of 
these equations are satisfying the constraints on quantum mechanical diffusion coefficients. 
Analytical expressions for the first two moments of coordinate and momentum are obtained 
by using the characteristic function of the Lindblad master equation. The master equa- 
tion is transformed into Fokker-Planck equations for quasiprobability distributions and a 
comparative study is made for the Glauber P representation, the antinormal ordering Q 
representation and the Wigner W representation. The density matrix is represented via 
a generating function, which is obtained by solving a time-dependent linear partial differ- 
ential equation derived from the master equation. Illustrative examples for specific initial 
conditions of the density matrix are provided. The solution of the master equation in the 
Weyl-Wigner-Moyal representation is of Gaussian type if the initial form of the Wigner 
function is taken to be a Gaussian corresponding (for example) to a coherent wavefunction. 
The damped harmonic oscillator is applied for the description of the charge equilibration 
mode observed in deep inelastic reactions. For a system consisting of two harmonic oscil- 
lators the time dependence of expectation values, Wigner function and Weyl operator are 
obtained and discussed. In addition models for the damping of the angular momentum are 
studied. Using this theory to the quantum tunneling through the nuclear barrier, besides 
Gamow's transitions with energy conservation, additional transitions with energy loss, are 
found. The tunneling spectrum is obtained as a function of the barrier characteristics. 
When this theory is used to the resonant atom-field interaction, new optical equations 
describing the coupling through the environment of the atomic observables are obtained. 
With these equations, some characteristics of the laser radiation absorption spectrum and 
optical bistability are described. 



1. Introduction 

In the last two decades, more and more interest arose about the problem of dissipation 
in quantum mechanics, i.e. the consistent description of open quantum systems [1-7]. The 
quantum description of dissipation is important in many different areas of physics. In 
quantum optics we mention the quantum theory of lasers and photon detection. There 
are some directions in the theory of atomic nucleus in which dissipative processes play 
a basic role. In this sense we mention the nuclear fission, giant resonances and deep 
inelastic collisions of heavy ions. Dissipative processes often occur also in many body or 
field-theoretical systems. 

The irreversible, dissipative behaviour of the vast majority of physical phenomena 
comes into an evident contradiction with the reversible nature of our basic models. The 
very restrictive principles of conservative and isolated systems are unable to deal with more 
complicated situations which are determined by the features of open systems. 

The fundamental quantum dynamical laws are of the reversible type. The dynamics of 
a closed system is governed by the Hamiltonian that represents its total energy and which 
is a constant of motion. In this way the paradox of irreversibility arises: the reversibility 
of microscopic dynamics contrasting with the irreversibility of the macroscopic behaviour 
we are trying to deduce from it. 

One way to solve this paradox of irreversibility is to use models to which Hamiltonian 
dynamics and Liouville's theorem do not apply, but the irreversible behaviour is clearly 
present even in the microscopic dynamical description. The reason for replacing Hamilto- 
nian dynamics and Liouville's theorem is that no system is truly isolated, being subject 
to uncontrollable random influences from outside. For this reason these models are called 
open systems. There are two ways of treating quantitatively their interaction with the out- 
side. One way is to introduce specific stochastic assumptions to simulate this interaction, 
the other one is to treat the open systems with the usual laws of dynamics, by regarding 
them as subsystems of larger systems which are closed (i.e. which obey the usual laws 
of dynamics with a well-defined Hamiltonian). The dissipation arises in general from the 
subsystem interactions with a larger system, often referred to as the reservoir or bath. 
The first of these two approaches has been used for the study of steady state transport 
processes in systems obeying classical mechanics. The second of the two approaches has 
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been mainly used in quantum mechanics. The main general result [2,8-10] is that under 
certain conditions the time evolution of an open system can be described by a dynamical 
semigroup $t(t > 0). For a closed finite system with Liouville operator the evolution 
operator is not restricted to nonnegative t. The importance of the dynamical semigroup 
concept is that it generalizes the evolution operator to open systems, for which there is no 
proper Liouville operator and no $t for negative t. The mathematical theory of dynamical 
semigroups has been developed in [2,11-15]. 

In Sect. 2 the notion of the quantum dynamical semigroup is defined using the concept 
of a completely positive map [14]. The Lindblad formalism replaces the dynamical group 
uniquely determined by its generator, which is the Hamiltonian operator of the system, by 
the completely positive dynamical semigroup with bounded generators. Then the general 
form of Markovian quantum mechanical master equation is obtained. 

The quantum mechanics of the one-dimensional damped harmonic oscillator represents 
a fundamental theoretical problem with applications in different domains of quantum op- 
tics, solid state physics, molecular and nuclear physics. In the present review paper the 
quantum harmonic oscillator is treated in the Lindblad axiomatic formalism of quantum 
dynamical semigroups. In Sect. 3 we give the fundamental constraints on quantum me- 
chanical diffusion coefficients which appear in the corresponding master equations [16,17]. 
On the basis of different representations it is shown that various master equations for 
the damped quantum oscillator used in the literature for the description of the damped 
collective modes in deep inelastic collisions or in quantum optics are particular cases of 
the Lindblad equation and that not all of these equations in literature are satisfying the 
constraints on quantum mechanical diffusion coefficients. In 3.1, by using the Heisenberg 
representation, explicit expressions of the mean values and variances are given [17]. In 
3.2 we solve the master equation with the characteristic function [18]. This function is 
found as a solution of a corresponding partial differential equation. By this method one 
can derive explicit formulae for the centroids and variances and, in general, for moments 
of any order. In 3.3 we explore the applicability of quasiprobability distributions to the 
Lindblad theory [19]. The methods of quasiprobabilities have provided technical tools of 
great power for the statistical description of microscopic systems formulated in terms of 
the density operator. In the present review the master equation of the one-dimensional 
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damped harmonic oscillator is transformed into Fokker-Planck equations for the Glauber 
P, antinormal Q and Wigner W quasiprobability distributions associated with the density- 
operator. A comparative study of these representations is made. The resulting equations 
are solved by standard methods and observables directly calculated as correlations of these 
distribution functions. We solve also the Fokker-Planck equations for the steady state and 
show that variances found from the P, Q and W distributions are the same [19]. In 3.4 we 
study the time evolution of the density matrix that follows from the master equation of the 
damped harmonic oscillator [20] . We calculate the physically relevant solutions (including 
both diagonal and off-diagonal matrix elements) of the master equation by applying the 
method of generating function. This means that we represent the density matrix with a 
generating function which is the solution of a time-dependent partial differential equation 
of second order, derived from the master equation of the damped harmonic oscillator. 
We discuss stationary solutions of the generating function and derive the Bose-Einstein 
density matrix as example. Then, formulas for the time evolution of the density matrix 
are presented and illustrative examples for specific initial conditions provided. The same 
method of generating function was already used by Jang [21] who studied the damping 
of a collective degree of freedom coupled to a bosonic reservoir at finite temperature with 
a second order RPA master equation in the collective subspace. In 3.5 the master equa- 
tion is given in the Weyl-Wigner-Moyal representation [17]. We show that the solution of 
the Lindblad equation in this representation is of Gaussian type if the initial form of the 
Wigner function is taken to be a Gaussian corresponding to a coherent wavefunction. 

In Sect. 4 we give some applications to nuclear equilibration processes. In 4.1 the 
charge equilibration mode is treated by a damped harmonic oscillator in the framework 
of the Lindblad theory [17]. It is shown that the centroids and variances of the charge 
equilibration mode observed in deep inelastic reactions can be well described by the cor- 
responding overdamped solutions [18,22]. In 4.2 we treat the damping of the proton and 
neutron asymmetry degrees of freedom with the method of Lindblad by studying the damp- 
ing of two coupled oscillators [23] . We present the equation of motion of the open quantum 
system of two oscillators in the Heisenberg picture and we derive the time dependence of 
the expectation values of the coordinates and momenta and their variances. We discuss 
also the connection with the Wigner function and Weyl operator. We demonstrate the 
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time dependence of the various quantities for a simplified version of the model, where the 
decay constants can be calculated analytically. 

In Sect. 5 simple models for the damping of the angular momentum are studied in the 
framework of the Lindblad theory of open systems [24,25]. We assume that the system is 
opened by the generators of the proper Lorentz group and the Heisenberg group. 

In Sect. 6 we study the quantum tunneling [26]. Considering the tunneling operator 
between the localized states, we solve Lindblad master equation in the second order ap- 
proximation of the theory of perturbations. We find that diffusion leads to additional 
tunneling processes around the Q-value, while friction determines only transitions at lower 
energies. Generally, this result is contrary to the majority of the previous studies [27,28] 
where, usually, it was concluded that the dissipation should decrease the tunneling rate. 

In Sect. 7 we study the resonant atom-field interaction [29]. From Lindblad's theorem, 
instead of the conventional optical Bloch equations, where only decay processes appear, 
we obtaine a more general system of equations where the atomic observables are coupled 
through the environment. Based on the new equations, experimental characteristics of the 
electromagnetic field absorbtion spectrum and of the optical bistability are described. At 
the same time, we predict that in some conditions, an energy transfer from the dissipative 
environment to the coherent electromagnetic field appears. 

The concluding summary is given in Sect. 8. 

2. Lindblad theory for open quantum systems 

The standard quantum mechanics is Hamiltonian. The time evolution of a closed 
physical system is given by a dynamical group Ut which is uniquely determined by its 
generator , which is the Hamiltonian operator of the system. The action of the dynamical 
group Ut on any density matrix p from the set ViTi.) of all density matrices of the quantum 
system, whose corresponding Hilbert space is H, is defined by 

for all t G ( — 00,00). We remind that, according to von Neumann, density operators 
p e ViTi) are trace class (Trp < 00), self-adjoint (p""" = p), positive (p > 0) operators with 
Trp = 1. All these properties are conserved by the time evolution defined by Ut- 

In the case of open quantum systems the main difficulty consists in finding such 
time evolutions for density operators p{t) = $t(p) which preserve the von Neumann 
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conditions for all times. From the last requirement it follows that must have the 
following properties: 

(i) $t(Aipi + A2P2) = Ai$t(pi) + A2$f(p2); Ai, A2 > with Ai + A2 = 1, 
(n) = $,(p)+, 

{in) > 0, 

(iv) Tr^tip) = 1. 

But these conditions are not restrictive enough in order to give a complete description 
of the mappings as in the case of the time evolutions Ut for closed systems. Even 
in this last case one has to impose other restrictions to Ut, namely, it must be a group 
Ut+s = UtUg. Also, it is evident that in this case Uo{p) = p and Ut{p) p m the trace 
norm when t — 0. For the dual group Uf acting on the observables A e B{7i), i.e. on the 
bounded operators on 7i, we have 

Ut{A) = e^^'Ae-^"'. 

Then Ut{AB) = Ut{A)Ut{B) and Ut{I) = /, where / denotes the identity operator on H. 
Also, Ut{A) A ultraweakly when t — > and Ut is an ultraweakly continuous mapping 
[2,9,13,14,30]. These mappings have a strong positivity property called complete positivity: 

^ B+Ut{A+Aj)Bj > 0, Ai, Bi e B{H). 

Because the detailed physically plausible conditions on the systems, which correspond to 
these properties are not known, it is much more convenient to adopt an axiomatic point 
of view which is based mainly on the simplicity and the succes of physical applications. 
Accordingly [2,9,13,14,30] it is convenient to suppose that the time evolutions $t for open 
systems are not very different from the time evolutions for closed systems. The simplest 
dynamics $t which introduces a preferred direction in time, characteristic for dissipative 
processes, is that in which the group condition is replaced by the semigroup condition 
[9-12,31] 

^t+S = ^t^3, t,S>0. 

The duality condition 

Tv{^t{p)A) = Tv{p^t{A)) (2.1) 
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defines the dual of acting on B{TC). Then the conditions 



and 

^t{I)=I (2.2) 

are equivalent. Also the conditions 

^ A (2.3) 

ultraweakly when t — > and 

^t{p) P 

in the trace norm when t — > 0, are equivalent. For the semigroups with the properties 
(2.2), (2.3) and 

A>0^ $t(A) > 0, 

it is well known that there exists a (generally bounded) mapping L, the generator of 

is uniquely determined by L. The dual generator of the dual semigroup $t is denoted 
by L: 

Tr{L{p)A) = Tr{pL{A)). 

The evolution equations by which L and L determine uniquely and respectively, are 
given in the Schrodinger and Heisenberg picture as 

!^ = i(«.W) (2.4) 

and _ 

'^=L(i,(A)). (2.5) 

These equations replace in the case of open systems the von Neumann-Liouville equations 

^ = -i[.,.KP)i 

and _ 
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For any applications Eqs. (2.4) and (2.5) are only useful if the detailed structure of the 
generator L{L) is known and can be related to the concrete properties of the open systems 
described by such equations. 

Such a structural theorem was obtained by Lindblad [14] for the class of dynamical 
semigroups $t which arc completely positive and norm continuous. For such semigroups 
the generator L is bounded. In many applications the generator is unbounded. 

A bounded mapping L : B{n) B{n) which satisfies L{I) = 0, L{A+) = L{A)+ and 

L{A+A) - L{A+)A - A+L{A) > 

is called dissipative. The 2-positivity property of the completely positive mapping $t: 

^t{A+A) > $,(A+)$,(yl), (2.6) 

with equality at t = 0, implies that L is dissipative. Lindblad [14] has shown that con- 
versely, the dissipativity of L implies that $t is 2-positive. L is called completely dissipative 
if all trivial extensions of L are dissipative. Lindblad has also shown that there exists a 
one-to-one correspondence between the completely positive norm continuous semigroups 
and completely dissipative generators L. The structural theorem of Lindblad gives the 
most general form of a completely dissipative mapping L [14]: 

Theorem: L is completely dissipative and ultraweakly continuous if and only if it is of 
the form 

L{A) = '-[H,A] + ^J2^V+[A,V,] + [V;,A]V,), (2.7) 

j 

where Vj, E,- V+Vj G B{n), H e B{n\,^, 

The dual generator on the state space (Schrodinger picture) is of the form 

L{P) = -it^' P] + 4 E([^^-P' + t^^' P^^])- (2.8) 
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Eqs. (2.4) and (2.8) give an explicit form for the most general time-homogeneous quantum 
mechanical Markovian master equation with a bounded Liouville operator: 

= L($,(p)) = [ff, $,(p)] + ^ X]([V;$t(p), ^+] + [VS, Up)V^]). 
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Talkner [32] has shown that the assumption of a semigroup dynamics is only apphcable in 
the hmit of weak couphng of the subsystem with its environment, i.e. for long relaxation 
times. 

We should like to mention that all Markovian master equations found in the literature 
are of this form after some rearrangement of terms, even for unbounded generators. 

It is also an empirical fact for many physically interesting situations that the time 
evolutions drive the system towards a unique final state p(oo) = limt_>oo ^t(p(0)) for 
allp(O) e v(n). 

From the 2-positivity property (2.6) of $t it follows that 

3 3 

or, by duality [17], 

T^{lt{p)Y,V;V^) > J2T:TiMp)V*)TriMp)Vj). (2.9) 

3 3 

The evolution equations of Lindblad are operator equations. The problem of finding their 
solutions is in general rather difficult. In cases when the equations are exactly solvable, 
these solutions give complete information about the studied problem and permit the cal- 
culation of expectation values of the observables at any moment. 

3. The one-dimensional damped quantum harmonic oscillator 
In this Section the case of the damped quantum harmonic oscillator is considered in 
the spirit of the ideas presented in the previous Section. The basic assumption is that the 
general form (2.8) of a bounded mapping L given by the Lindblad theorem [14] is also 
valid for an unbounded completely dissipative mapping L: 

L{P) = ~[H,p] + ^Y.([VjP,V;] + [V„pV+]). (3.1) 

3 

This assumption gives one of the simplest way to construct an appropriate model for a 
quantum dissipative system. Another simple condition imposed to the operators H, Vj, 
is that they are functions of the basic observables q and p of the one-dimensional quantum 
mechanical system (with [g',p] = %hl^ where / is the identity operator on 'H) of such kind 
that the obtained model is exactly solvable. A precise version for this last condition is 
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that linear spaces spanned by the first degree (respectively second degree) noncommutative 
polynomials in p and q are invariant under the action of the completely dissipative mapping 
L. This condition implies [16] that Vj are at most first degree polynomials in p and q and 
H is at most a second degree polynomial in p and q. 

Because in the linear space of the first degree polynomials in p and q the operators p 
and q give a basis, there exist only two C-linear independent operators Vi , V2 which can 
be written in the form 

Vi = ttip + biq, i = l,2, 

with ai,bi — 1,2 complex numbers [16]. The constant term is omitted because its contri- 
bution to the generator L is equivalent to terms in H linear in p and q which for simplicity 
are assumed to be zero. Then H is chosen of the form 

H = Ho + f^{pq + qp), H, = ^p'+'^q\ (3.2) 

With these choices the Markovian master equation can be written: 

dp % % % 

= -^[Ho,p] - ^{>^ + IJ')[q,PP + Pp] + ^{>^ - IJ')\P, PQ + QP] 

hp]] - ^b, b,p]] + ^{[1, \p,p]] + \p, hp]])- (3.3) 

n n n 

Here we used the notations: 

i=l,2 3 = 1,2 i=l,2 j=l,2 

where Dpp,Dqq and Dpg are the diffusion coefficients and A the friction constant. They 
satisfy the following fundamental constraints [17]: 

i) Dpp > 0, 

ii) Dqq > 0, (3.4) 

in) DppDqg - Dpq^ > 

Indeed, inequalities i) and ii) follow directly from the definitions of the coefficients and Hi) 
from the Schwartz inequality. The equality p, = X is a, necessary and sufficient condition 
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for L to be translatfon invariant [16]: [p, L(p)] = L([p, p]). In the foUowing general values 
for A and ji will be considered as in [16]. 

Introducing the annihilation and creation operators 

1 i 1 i 

a= ^=={^/mwq+ —=p), a+ = -^{y/mwq ^=p), (3.5) 

V2?i vrmo ^/2h \Jmw 

obeying the commutation relation [a, = 1, we have 

i^o = ^a;(a+a+ -) (3.6) 
2 



and the master equation has the form 

^ - \Dx - ii){pa+a+ - a+pa+) + -(Di + ii){a+a+p - a+pa"^) 



where 



dt 

1 1 
+ -{D2 — A — iuj){a^ pa — paa^) + -{D2 + A + iio){apa'^ — a^ap) + /i.e., (3-7) 

Di = \{mujDqq-^ + 2iDpq), D2 ^ UmuDgg + (3.8) 



In the literature, master equations of the type (3.3) or (3.7) are encountered in models 
for the description of different physical phenomena in quantum optics, in treatements of 
the damping of collective modes in deep inelastic collisions of heavy ions or in the quantum 
mechanical description of the dissipation for the one-dimensional harmonic oscillator. In 
the following we show that some of these master equations are particular cases of the 
Lindblad equation. 

1) The master equation of Dekker for the damped quantum harmonic oscillator [5,33- 
36] supplemented with the fundamental constraints (3.4) obtained in [36] from the con- 
dition that the time evolution of this master equation does not violate the generalized 
uncertainty principle [37] at any time, is a particular case of the Lindblad master equation 
(3.7) when /i = A. If the constraints (3.4) are not fulfilled in the case p, = X, then the 
uncertainty principle is violated. 

2) The quantum master equation considered in [38,39] by Hofmann et al. for treating 
the charge equilibration process as a collective high frequency mode is a particular case of 
the Lindblad master equation (3.3) if A = 7(a;)/2m = p, Dqq = Dpq = 0, Dpp = 7(a;)T*(a;), 
but the fundamental constraints (3.4) are not satisfied. 
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3) For the quantum master equation considered in [40] for the description of heavy ion 
collisions we have \ = ji = 7/2, Dpp = D, Dqq = 0, Dpg = Dqp = —d/2 and consequently 
the fundamental constraints are not fulfilled. 

4) In [41], Spina and Weidenmiiller considered two kinds of master equations for 
describing the damping of collective modes in deep inelastic collisions of heavy ions. 

(i) The first equation can be obtained from Eq. (3.3) by replacing Hq by 

Hq — Amojq'^/2 + f{t)q and setting \ = jj, = r/2, Dpp = D/2, Dqq = and Dpq = 
Dqp = B/2. Then the constraints (3.4) are not satisfied. 

(ii) The second master equation is obtained from (3.3) by setting Hq — AmLijq'^/2 
—Ap^/{2mu) + f{t)q instead of Hq and = Tp = A,// = 0,Dpp = Dp/2, Dqq = 
Dfi/2, Dpq = 0. Then last condition (3.4) is satisfied for all values of the parameters. 

5) The master equation for the density operator of the electromagnetic field mode 
coupled to a squeezed bath [42,43] can be obtained from the master equation (3.7) if we 
set /X = 0, A = 7, Di = 27M, D2 = 7(27V + 1) (in [42]: 7/2 instead of 7 and N instead of 
AT). 

6) The master equation for the density operator of a harmonic oscillator coupled to 
an environment of harmonic oscillators considered in [44-47] is a particular case of the 
master equation (3.7) if we set A = /U = 7, Dqq = Dpq — 0, Dpp — 27(77, + l/2)mwTi and the 
fundamental constraints (3.4) are not fulfilled. 

7) The master equation written in [48] for different models of correlated- emission 
lasers can also be obtained from the master equation (3.7) by setting 

Di+ iJL = 2A4, Di- fjL^ 2A3, D2 + \ + ioJ^ 2A2, D2- X-iuo ^ 2Ai. 

8) Two master equations were introduced by Jang and Yannouleas in [21,49], where 
the nuclear dissipative pocess is described as the damping of a collective degree of freedom 
coupled to a bosonic reservoir at finite temperature. The resulting RPA master equation 
within the observed collective subspace was derived in a purely dynamical way. 

(i) The master equation written in [49] in the resonant approximation (rotating-wave 
approximation) can be obtained as a particular case of the Lindblad master equation (3.7). 
For this one has to set 

Dpp = m^oj^Dqq, Dpq = ii = 0, = (2 < n > +l)r, A = ^, 
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where < n > is the average number of the RPA collective phonons at thermal equilibrium 
and r the decay constant. The fundamental constraints (3.4) are fulfilled in this case. 

(ii)The master equation derived in [21] in order to extend the calculations carried out 
in [49] with the before-mentioned master equation, can also be obtained as a particular 
case of the master equation (3.7) by setting 

Dqq = Dpq = 0, Dpp = ^!^(2 < n > +l)r, i_i = X = ^ 

or D2 = —Di = (2 < n > -|-l)r/2. In this case the fundamental constraints (3.4) are not 
fulfilled. 

As a conclusion, we can say that on various occasions the master equation for X = /i 
has been tested in the sense that it does not violate the uncertainty principle in the final 
steady state, but a possible transient violation of the uncertainty principle has been noted 
[5,36]. The above mentioned models can be divided into three types. The first type are 
models which fulfill the constraints (3.4). The corresponding master equations belong 
to the cathegory of equations described by the Lindblad theory. These models obey the 
uncertainty principle. The second type are models with X — jj. which do not fulfill the 
constraints (3.4) (for example model 2). According to [36] these models show a dynamical 
violation of the uncertainty principle. The third type are models with X^ jj. which also do 
not fulfill the constraints (3.4). These models and those of the second type do not belong 
to the cathegory of master equations treated in the Lindblad theory. For the models of the 
third type we can not prove that they do not conserve the uncertainty principle in time. 

3.1 The master equation in the Heisenberg picture 

The following notations will be used: 

aq{t)^Tr{p{t)q), 

ap{t) = Tv{p{t)p), 
aqq = T:T{p{t)q^)-al{t), 
app = Tr{p{t)p'')-al{t), 
apq{t) = T,{p{tf-^±^)-ap{t)aq{t) 
12 



(3.9) 



In the Heisenberg picture the master equation has the following symmetric form: 

+ - ^){QmA),p] + mA),p]q) - ^[q, [q, MA)]] 

^''■\p, \p, MA)]] + [g, MA)]] + [g, [p, MA)]]). (3.10) 



Denoting by ^ a selfadjoint operator we have 

aA{t) = Tr {p{t) A), a AA{t) = TT{p{t)A^) - a\{t). 

It follows that 

= Tr{L{p{t))A) = T^iPmiA)) (3.11) 

and 

= Tr{L{p{t)A') - 2^^aA{t) = Tr{p{t)L{A')) - 2aA{t)Tr{{p{t)L{A)). (3.12) 

An important consequence of the precise version of solvability condition formulated at 
the beginning of the present Section is the fact that when A is put equal to p or g in 
(3.11) and (3.12), then dap{t)/dt and daq{t)/dt are functions only of cTpit) and (Jq{t) and 
dapp{t) /dt, daqq{t) /dt and dapq{t)/dt are functions only of app{t),aqq{t) and apq{t). This 
fact allows an immediate determination of crp{t), aq{t) , app{t) , aqq{t) and crpq{t) as functions 
of time. The results are the following: 

^ = -(X-^^)aq{t) + ^apit), 
dt m 

'^''P^'^ = -moo'^aqit) - (A + p)ap{t) (3.13) 



and 



dt 



= -2(A - p)aqq{t) + '^(Tpqit) + 2Dqq, 

d(T 

= -2(A + p)app{t) - 2muj^apq{t) + 2Dpp, (3.14) 

= -muj''aqq{t) + —(Tpp{t) - 2Xapq{t) + 2Dpq. 
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All equations considered in various papers in connection with damping of collective modes 
in deep inelastic collisions are obtained as particular cases of Eqs. (3.13) and (3.14), as we 
already mentioned before. 

The integration of Eqs. (3.13) is straightforward. There are two cases: aju > u> 
(overdamped) and b) n < to (underdamped). 

In the case a) with the notation u"^ = /j,"^ — a;^, we obtain [17]: 

ag{t) = e~^*((coshzyt + - siBhut)aq{0) + — siiahi'tap{0)), 

2 

ap{t) = e~^*(-^^ sinhzyta-q(O) + (coshz^t - ^ smhut)ap{0)). (3.15) 

If X> v, then crq(oo) = (Tp(oo) = 0. If A < i^, then crg(oo) = (Jp(oo) = oo. 
In the case b) with the notation fl^ = uj'^ — /j,^ , we obtain [17]: 

aq{t) = e~^\{cosnt+ ^sinrit)crq(O) + smntap{0)), 

O 

ap{t) = e-^*( —smntaq{0) + (cosOt - ^smnt)ap{0)) (3.16) 

and CTq (oo) = CTp (oo) = 0. In order to integrate Eqs. (3.14), it is convenient to consider the 
vector 

/ mU)(Tqq{t) 

X{t) = I app{t)/muj 

Then the system of equations (3.14) can be written in the form 

dX{t) 
dt 

where R is the following 3x3 matrix 



RX{t) + D, 



-2(A-^) 2u; 

R = \ -2(A + /i) -2io 

—oj (jj — 2A 



and D the following vector 



2mu)Dqq 



D = I 2Dpp/muj 



2D 

P9 
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Then there exists a matrix T with the property T'^ = I where / is the identity matrix and 
a diagonal matrix K such that R = TKT. From this it follows that 



X{t) = (Te^*T)X(0) + T(e'^* - I)K-'TD. (3.17) 

An interesting observation is that the time dependence of the variances aqq{t), app{t), 
(Tpq{t) decomposes in a classical part given by Te^*TX(0) and a quantum part given by 
T(e^* - I)K-^TD. Exactly this quantum part governs the asymptotic behaviour of the 
variances when t — > oo [17]. 

In the overdamped case {fj, > oj) the matrices T and K are given by 



1 



and 



T — — \ (1 — V iJL-\- V 2u) 
—ijj —ijj 

-2(A -v) 

K=\ -2{\ + v) 

-2A 

with = - J^. 

In the underdamped case {\i < to) the matrices T and K are given by 

I / /i + ifl li — iVt 2u) 

T = — — \ n-ifl n + in 2u 

2ziM r, 

\ —U3 —UJ —IjJL 

and 

(-2{X-iVt) 
K=\ -2{\ + in) 

\ -2A 

with Q."^ = oj^ - fi^ . From (3.17) it follows that 

X{oo) = -(TK-^T)D = -R-^D (3.18) 

(in the overdamped case the restriction X > u is necessary). Then Eq.(3.17) can be written 
in the form 

X(t) = {Te^*T)(X(0) - X(oo)) + X(oo). (3.19) 



Also 

dX{t) ^ ^TXe^'T){X{Q) - X{oo)) = R{X{t) - X{oo)). 



dt 
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The formula (3.18) is remarcable because it gives a very simple connection between the 
asymptotic values (t — > oo) of aqq{t), app{t), crpq{t) and the diffusion coefficients Dgg, Dpp, 
Dpq. As an immediate consequence of (3.18) this connection is the same for both cases, 
underdamped and overdamped, and has the following explicit form: 



1 

2(ma;)2A(A2+a;2 -//2)' 



^^^M = nf^,.^2^f^2 , ,,2 , ((ma;)'(2A(A + /x) + a;^)I^,g 



+uP'Dpp + 2mu;^(A + ii)Dpq), 
(Tpp{oo) = 2A(A2 +1,2 _ ^2) ((^) V£)gg+(2A(A-//) W)i^pp-2W(A-/x)£)pg), (3.20) 

(7pq(00) = 2y^A(A2 + - ^2) (~(^ + lAi'^T^YDqq + (A - //)i:'pp + 2m(A^ - ll'^)Dpq). 

These relations show that the asymptotic values a"qq((X)), a"pp(oo), a"pq(oo) do not depend 
on the initial values (Tqq(O), (Tpp(O), (Tpq(O). In other words, 

_^ / 2A(A + //) + a;^ 2a;(A + //) 

4A(A2+a;2-^2) -(A + ^)a; (A -^o. 2(A^ - /) 

Conversely, if the relations D = —RX{po) are considered, then 

Dqq = (A - lj)aqq{00) - — CTpq(oo), 

-Dpp = (A + ij)app{oo) + ■muj^apq{oo), (3-21) 

-Dpq = ^(??^a;^(7qq(oo) - — Crpp(oo) + 2Acrpq(oo)). 

Hence, from (3.4) fundamental constraints on a"gq(oo), a"pp(oo) and (7pq(po) follow: 

Dqq = (A - IJi)aqq{00) - — C7pg(oo) > 0, 

Dpp = (A + /x)(jpp(oo) + muj'^apq{oo) > 0, 
DqqDpp - Dig = (A^ - //^)a-qq(oo)c7pp(oo) - cu^cTpg (oo) + (A - ii)mu)'^ aqq{oo)apq{oo) 
-^^^^^J^(^Pp{oo)apq{oo)-^{mu;^falq{oo)-^^ 

A A^n^ 

—muj^ Xaqq{oo)apq{oo) H (Tpp(oo)(Tpq(oo) > — — . (3.22) 
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The constraint (3.22) can be put in a more clear form: 



4(A^ + a;^ - //^)(crgq(oo)crpp(oo) - crpq(oo)^) 

-(ma;^a-qq(oo) + — (jpp(oo) + 2/X(Tpg(oo)) > %^\^. (3.23) 

If < a; (the underdampcd case), then + a;^ — //^ > A^. If // > w (the overdamped 
case), then < A^ + a;^ — //^ < A^ (A > iv) and the constraint (3.23) is more strong than 
the uncertainty inequahty aqq{oo)app{oo) — (7^^(00) > Ti? /A. By using the fact that the 
hnear positive mapping B{T-L) — > C defined by ^ — > Tr(p^) is completely positive (hence 
2-positive), in [17] the following inequality was obtained: 

Dqqapp{t) + Dppaqq{t) - 2Dpqapq{t) > 

We found this inequality, which must be valid for all values of t e (0, 00), as the sufficient 
condition that the generalized uncertainty inequality [37] 



4 



is fulfilled at any moment if the initial values crgg(O), crpp(O) and (Jpq{Q) for t = satisfy 
this inequality. A restriction connecting the initial values (Tqq(O), (Tpp(O), (7pq(0) with the 
asymptotic values crpp(oo), crqq(oo), crpq(oo) is also obtained if the values Dqq, Dpp and Dpq 
are expressed by Eq. (3.21) in terms of (7pp(oo), (7qq(oo), (7pq(oo) : 

Dqq(7pp{0) + Dppaqq{0) - 2Dpqapq{0) > 

More explicitly 

A(crgq(oo)crpp(0) + (Tpp{oo)aqq{Q) - 2(7^5 (oo)crpg(O)) - //(crgq(oo)crpp(O) - Crpp(oo)(7gq(0)) 
1 

((^Pg(OO)o-pp(0)-(Tpp(OO)(Tpg(0))+ma;^((Tpg(OO)(Tgg(0)-(Tgg(OO)(Tpg(0)) > — " . (3.24) 

If the asymptotic state is a Gibbs state 

PG(oo) = e-^°/^^/IV(e-^°/'=^), 
17 



then 
and 

„ A + u. ^ hu> ^ X — a h ^ hiv ^ ^ 

Dr,r, = — - — nmucota——, Daa = — 7: cotn— — , Dj,o = (3.26) 

and the fundamental constraints (3.4) are satisfied only if A > /x and [16]: 

(A^-/^^)(coth|^f >A^ 
If the initial state is the ground state of the harmonic oscillator, then 

Then (3.24) becomes 

A((Taq(oo)?7iu; H — — ) — //((Tqq(oo)ma; — ) > nX. 

mw rriLO 

For example, in the case (3.25), this implies coth(7ic<;/2/cT) > 1 which is always valid. 

Now, the explicit time dependence of aqq{t), (Tpp{t), (Jpqit) will be given for both under- 
and overdamped cases. From Eq. (3.19) it follows that in order to obtain the explicit time 
dependence it is necessary to obtain the matrix elements of Te^^T. In the overdamped 
case (// > u;, z/^ = — u'^) we have 



Te^^T = '^^^^ I a2i 022 0^23 

032 033 



e 



2Xt / CtU 0^12 0^13 




1 



with 

O'li = (a*^ + ^'^) cosh 21/t + 2/11' sinh 21/t — a;^, 

0-12 = (a*^ — i^^) cosh2i>'t — a;^, 

013 = 2c(;(/iCOsh2z/i + ^'sinh2z>'t — //), 

^21 = (a*^ — '^^) cosh2i/t — w^, 

022 = (At^ + z^^) cosh 2z/t - 2niy sinh 2zyt - w^, (3.27) 

023 = 2a;(//cosh2z^t — ^'sinh2i/t — //), 

031 = —uj {jJL cosh 2z>'t + V sinh 2z>'t — //) , 

032 = —Lo{iJL cosh.2vt — vsmh.2vt — /i), 
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a33 = -2(a;^cosh2i/t - //^). 

In the underdamped case {/i < a;, = o;^ — //^) we have 

-2At / hi bi2 feisX 
\ 031 032 033 / 

with 

611 = (^2 _ 20t - 2;uO sin 20^ - u;2, 

h2 = (/i^ + ^^^) cos20t - a;^ 

613 = 2a;(/icos2f]t — f]sin2f]t — //), 

b2i = (/i^ + fi^) cos 2nt - J^, 

b22 = iiJ,'^ -O.'^) cos 2nt + 2fj,nsm2nt- uj^, (3.28) 
623 = 2a; cos 2Qt + sin 2^},t — //), 

631 = —a; (/X cos 2nt — Q sin 2Qt — /u) , 

632 = —a; {fi cos 2f]t + sin 2f]t — //) , 

633 = -2{uj^cos2nt - fj,'^). 

As was shown by Lindblad in [16] the equations of motion (3.10) written for Weyl 
operators (A = VF(^, ry) = e^*/^^*^''^"^^)) can be integrated in a very simple and elegant way. 
This fact is important from both points of view: practical and theoretical. From practical 
point of view, it gives in a new way explicit formulas for agg(t), (7pp{t), 0'pq{t) and moreover, 
it gives the action of the dynamical semigroup $t generated by (3.10) on any polynomial in 
the noncommutative variables p and q. Because ^t{AB) ^ $t(yl)$t(S) it is not sufficient 
to know ^t(p) and ^t(q') as in the case of dynamical groups Ut{AB) = Ut{A)Ut{B), where 
evidently it is suflBcient to know the action of Ut on p and q in order to know this action 
on all noncommutative polynomials. From theoretical point of view, the explicit action 
of $t on Weyl operators W{^,ri) allows, as was shown by Lindblad in the Appendix of 
[16], directly to prove that the semigroup $t is indeed a semigroup of completely positive 
mappings. This assertion cannot be considered as a consequence of the structural theorem 
of Lindblad because L is unbounded. 

We can write as Lindblad [16], 
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and then the master equation 

d^t{W{C,rj)) 



L{^tiWiC,rj))) 



dt 

is equivalent with the following differential equations 

at m 



dt 



dgjt) ^ ^2u\ , n „2 

dt 



where ^(0) = ^, r]{Q) — rj and g{0) — 0. The explicit determination of ^t{W{$,,r])) as a 
function of t reduces to the integration of these differential equations. We obtain [17]: 

= ait)m + Pit)vio), vit) = iit)m + mm, 

where 

«(^) Pit) \ ^ g-At ( cosh lyt- ^ sinh lyt - ^ sinh lyt 
7(t) 5{t) J \ in^ sinh lyt cosh z^t + ^ sinh z^t 

From 

^t{W{^{0), r/(0))) = W{^{t), r/(t))e^W = e(Vft)('7Wg-^Wp)+3W^ (3.29) 
\m=v(0)=o - -^P 

and 

k(o)=,(o)=o - ^g, 

the action of $t on q and p, respectively, is obtained: 

Mq) = - /3{t)p, $t(p) = -7{t)q + a{t)p. 

As an immediate consequence of these relations we have [^t{Q)i^t{p)\ = ihle~'^^^*. From 
(3.29) and 

5^2(0) k(0)=r,(0)=0 - -^q 
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and 

d^w{m,vm, 1 



it follows, respectively, 
and 



where 



/ 2B{t)mu \ 
2A{t)/mw ^T{e^^ - I)K-^TD. 



muaqq{t) Cfpq{t) 
<^pq{t) (Tpp{t)/mU! 



Let us denote by a{t) the correlation matrix 

a{t)- 

Now by definition, for any state p E 'D{H) : 

aqq{t) =Tr{p^t{q^)) - {Tr{p^t{q))f 
appit) = TV(p$,(p2)) - (IV(p$,(p)))^ 



and 



Using the relations (3.20) and (3.21) the foUowing interesting formulas are obtained 

(A^ + a;2 - (7(00) = - det L> + -y^i^Dpp + -^D^q + pDpqf (3.30) 

and 

(A^ + o;^ - //^) det (j(oo) = - det D + (- — a-pp(oo) H — (Tqq(oo) + papq{oo))'^, (3.31) 

where 

\^ Dpp/muj 
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Comparison of (3.30) with (3.31) gives 



— (jpp(oo) + ^-(Tgg(oo) + iiapq[Qo) = ji^^pp + ^"^gg + I^Dpq). (3.32) 

But, the left-hand side of (3.32) is exactly the asymptotic mean value of the energy of the 
open harmonic oscillator. Hence (3.32) gives the value of E{oo) as a function of diflFusion 
coefficients, 

11 muj'^ 

E{00) = ji^Dpp + -^Dqq + llDpq). 

Another expression for E{oo) which also follows from (3.30) and (3.31) is 

E{oo) = ((A2 + a;^ - /x^) det a{oo) - ^ det bf''. 
The reality condition for E{oo) implies 

det (t(oo) > — — ^ 5 det D > 



4(A2 + a;2-^2) - 4(A2 + a;2 - /^s) 

which is more restrictive for the overdamped case (a; < //) than the generalized uncertainty 
inequality deta(oo) > 

3.2 The method of the characteristic function 

Instead of solving the master equation (3.7) directly, we first introduce the normally 
ordered quantum characteristic function x(-^) ■^*) t) defined in terms of the density operator 
P by 

X(A, A*, t) = Tr[p(t) exp(Aa+) exp(-A*a)], (3.33) 

where A is a complex variable and the trace is performed over the states of system. Sub- 
stituting Eq.(3.33) into the master equation (3.7) and using the operator relations 

aexp(Aa"'') = exp(Aa+)(a + A), a+exp(-A*a) = exp(-A*a)(a+ + A*), 

Tr[p(t) exp(Aa+) exp(-A*a)(a+ + A*)] = ^ax, Tr[p(t) exp(Aa+) exp(-A*a)a] = -^a^x 
or applying the rules: 

d d d d 

p ^ X, op ^ ~aA^^' ""^P ^^dA~ ^ ^~dA^ ^ ^""^ ^ dX^' 
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the following partial differential equation for x is found [18]: 

{dt + [(A - iu;)A + MI^A + [(A + iu;)A* + M]9a'}x(A, A*, t) 

= {L|A|2 + CA^ + C*A*2}x(A, A\t), (3.34) 

where 

L = X-D2, C=^{f, + D*). 

We consider the state of the system initially to be a superposition of coherent states. 
Coherent states |q; > of the harmonic oscillator are minimum uncertainty states having 
the mean coordinate < q > and mean momentum < p > given by 

< q >—< alqla >— \/ — ^Rect, < p >—< ablo: >= \/2hm(jjIma. (3.35) 
V mco 

Consequently, we take as the initial density operator 



p(0) = J dadpN{a,P)\a >< P\. 



The quantum characteristic function corresponding to the operator jo; >< P\ is given from 
Eq. (3.33) by 

X=</3|a >exp(A/3* - A*a). (3.36) 
We look for a solution of (3.34) having the exponential form 

x(A,A*,t) = J dad(3N{a,(3) < (3\a > exp[A{t)A + B{t)A'' + f{t)A^ + f*{t)A*^ + h{t)\A\^]. 

(3.37) 

The form of the solution (3.37) is suggested from the fact that the left-hand side of Eq. 
(3.34) contains first-order derivatives with respect to the time and variables A and A* and 
is symmetric with respect to complex conjugation. The functions A{t), B{t), f{t) and h{t) 
depend only on time. Corresponding to Eq. (3.36), these functions have to satisfy the 
initial conditions 

A(0) = S(0) = a, /(O) = 0, /i(0) = 0. (3.38) 

When we introduce the function (3.37) into Eq. (3.34) and equate the coefficients of equal 
powers of A and A*, we get the following two systems of linear differential equations of 
first order with constant coefficients: 

+ (A - iuj)A{t) + iJiB{t) = 
23 



dt 
and 



^^^^'^ + ^A{t) + {X + iuj)B{t) = (3.39a) 



+ 2XR{t) + 2ul{t) + nh{t) = ReC 
dl{t) 



dt 

dh{t) 



+ 2\I{t) - 2uR{t) = ImC (3.396) 
+ 4iiR{t) + 2Xh{t) = L, 



dt 

where R{t) = Re/(t),/(t) = Im/(t) with the initial conditions R{0) = 1(0) = h{0) = 0. 
Subject to the initial conditions (3.38), the homogeneous system (3.39a) has the solution 
[18]: 

A{t) = u{t)P* - v{t)a, B{t) = -u*{t)a + v{t)P*, (3.40) 

where 

u{t) = ^[exp(-/i_t) +exp(-/x+t) + — — — (exp(-^+t) - exp(-/x_t))], 

v{t) = (exp(-//_t) - exp(-/x+t)). (3.41) 

fj,- - fj.+ 

The eigenvalues are given by 



A ± - a;2, 7 = _ (3 42) 



The system (3.396) has the eigenvalues — 2A, — 2(A ± -y///^ — o;^) = — 2//± and in order to 
integrate it we apply the same method as for the system (3.14) in the preceding Subsection. 
We obtain: 

f{t) = ^ exp(-2//+t)(7 - - ^ exp(-2//_t)(7 + iu) - exp(-2At) + /(oo), 

/i(t) = Mexp(-2At) + Nexp{-2ii_t) + Pexp(-2//+t) + h{oo). (3.43) 
Here M, A^, P, /(oo) and /i(cxd) are constants given by [18]: 

M = ^i,T^C +f).N= ^^^(7R.C - .toC - f ), 
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and the asymptotic values connected with the diffusion coefficients -Dqq, -Dpp and Dpq are: 

2(AReC - c^ImC) - LfjL 



R{oo) = 



4(A2 - 72) 



'^^^ = ^W^) ' ^'-"'^ 

, , , L(A2 + 00^) - 2fi{XReC - ojIuiC) 
= 2A(A2-72) • 

By knowing the characteristic function (3.37), (3.40)-(3.43) corresponding to the initial 

density operator which represents a superposition of coherent states, it is easy to obtain 

explicit formulae for the moments: 

< a+^{t)a-{t) >= Tr[a+-^{t)a-{t)p{t)] = {-^r g^^^g^^ xi^, A\t)U=A*=o. 

In the following, we take the density operator p in the coherent state representation 

p(0) = J P{a)\a >< a\d^a, 

where -P(ct) is the diagonal or Glauber P distribution and d'^a = dKeadlma. The integra- 
tion covers the entire complex a plane. Then the characteristic function (3.37) becomes: 

X(A, A*, = J d^aP{a) exp[{ua* - va)A + {-u*a + va*)A*] expf/A^ + f*A*^ + h\Af]. 

Let us assume that the damped oscillator is at t = prepared in a pure coherent state, 
say \ao >, corresponding to P{a) — S(Rjea — Rjeao)6{Ima — Imcuo)- One has 

X^^\A, A*,t) = exp[(«a* - vao)A + {-u*ao + va*o)A*] exp[fA^ + f*A*^ + h\Af]. 

The first moments are given by 

< a^{t) >= — — — |a=a*=o = uaQ - vao, < ait) >= gA* 1^^^*^° = u ao - vaQ. 

Then, with the notations (3.9), using (3.35) and the transformations 
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for the displacement operator q{t) and the momentum operator p{t) of the oscillator, we 
obtain the following mean values: 



with u,v given by (3.41), (3.42). There are two cases: 

a) the overdamped case: /i > u , = /i'^ — uj^ , 'j = v; then 

u{t) = exp(— At)(coshz>'t + — sinhi/i), v{t) = — — exp(— At) sinhi/t (3.45) 

and Cqit), cfpit) take the previous form (3.15); 

h) the underdamped case: jji < ui^Q? = uj"^ — /j?, 7 = iQ; then 

u{t) = exp(-At)(coshr]t+ ^sinhm), v{t) = exp(-At) sinhl^t (3.46) 

and CTqit), (7p{t) take the previous form (3.16). 
For the variances one finds: 

< a\t) >= ^^^l^ \k=k*=o = («*«o - valf + 2/*, 

< a+\t) >= g^,^ ^ |a=a*=o = (««o - ^«o)' + 2/, (3.47) 

< a+(t)a(t) >= — gAM* 1^^^*^° ^ ^^"0 " vao){u*aQ - va^) - h. 

Then the relations (3.9) will give us the explicit time dependence of the variances (7qq{t), 
(Tpp(t), apqit). The asymptotic values of these variances are given by the following expre- 
sions [18]: 

Ti 1 

(7qq{00) = {f + f* - h + -)|t^oo, 

c7pp(oo) = -hmu{f + f* + h- ^)\t^oo, 

apq{oo) = ih{f - /*)|t_oo- 

With /(oo) = i?(oo) + iI{oo) and by using the formulas (3.44) for i?(oo), I{oo), h{oo), the 
asymptotic values of the variances take the same form (3.20) as in the preceding subsection, 
as expected. 
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With the relations (3.47), the expectation value of the energy operator can be calcu- 
lated: 



3.3 Quasiprobability distributions for the damped harmonic oscillator 

The methods of quasiprobabilities have provided technical tools of great power for 
the statistical description of microscopic systems formulated in terms of the density op- 
erator [50-53]. The first quasiprobability distribution was introduced by Wigner [54] in 
a quantum-mechanical context. In quantum optics the P representation introduced by 
Glauber [55,56] and Sudarshan [57] provided many practical applications of quasiprobabil- 
ities. The development of quantum-mechanical master equations was combined with the 
Glauber P representation to give a Fokker-Planck equation for the laser [58,59]. One use- 
ful way to study the consequences of the master equation for the one-dimensional damped 
harmonic oscillator is to transform it into equations for the c-number quasiprobability dis- 
tributions associated with the density operator. The resulting differential equations of the 
Fokker-Planck type for the distribution functions can be solved by standard methods and 
observables directly calculated as correlations of these distribution functions. However, 
the Fokker-Planck equations do not always have positive-definite diffusion coefficients. In 
this case one can treat the problem with the generalized P distribution [51]. 

First we present a short summary of the theory of quasiprobability distributions 
[53,60]. For the master equation (3.7) of the harmonic oscillator, physical observables 
can be derived from the expectation values of polynomials of the annihilation and creation 
operators. The expectation values are determined by using the quantum density operator 
p. Usually one expands the density operator with the aid of coherent states, defined as 
eigenstates of the annihilation operator: a\a >— a\a >. They are given in terms of the 
eigenstates of the harmonic oscillator as 



with the normalization | < P\a > P = exp(— ]q;— /3|^). In order to solve the master equation 
(3.7) we represent the density operator p by a distribution function over a c-number phase 
space. The chosen distribution function, introduced in [60], is defined as follows: 






(3.48) 
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with the characteristic function 

x{A,s) = Tr[pD{A,s)], 
where D{A, s) is the displacement operator 

^(A, s) = exp(Aa+ - A*a + ^s|Ap). 
The interval of integration in Eq. (3.48) is the whole complex A plane. Because of 

8'^ [a) = j exp{aA* - a*A)(fA, 

the distribution function $(q;,s) is the Fourier transform of the characteristic function. 
Since the density operator is normalized by Trp = 1, one obtains the normalization of 

In this paper we restrict ourselves to distribution functions with the parameters s = 1, and 
— 1. These distribution functions can be used to calculate expectation values of products 
of annihilation and creation operators. For that purpose we first expand the displacement 
operator in a power series of the operators a and a'^: 



oo oo 



DiA, .) = E E 1,., ^ K-^^n.- (3.49) 



m\n\ 

m=U n=(J 



The parentheses with the index s indicate the special representations of the polynomials 
depending on s. For example for n = m = 1 we have the s-ordered operators: 

{a~^a}s=i = a^a, {a~^a}s=Q = -{a'^a + aa'^), {a^a}s=-i = aa"*". (3.50) 

Expectation values of the s-ordered operators can be calculated as follows: 

< {a+-a"}. >= TV[p{a+-a"},] = ^[p(|^)"^(-^)^I^(A, s)|a=o] 

For the last step we apply the inverse relation to Eq. (3.48): 

x(A,s) = y $(Q;,s)exp(AQ;* - A*Q;)d^Q;. 
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In the following we discuss the distribution functions for s = 1,0 and —1 in more detail. 
For s = 1 we obtain the Glauber P function [56,61], for s = the Wigner function [54] 
and for s = — 1 the Q function [60]. For s = 1 we have 



Then the s ordering in Eq. (3.49) corresponds to normal ordering. Since the Glauber P 
function is the Fourier transform of the characteristic function 



it follows from Eq. (3.48) that the distribution $(«, 1) is identical to the P function. 
This function is used for an expansion of the density operator in diagonal coherent state 
projection operators [55-57,62]: 



Calculating the expectation values of normally ordered operator products we obtain the 
relation 



from which we again derive P{(y) = $(a, 1). Despite the formal similarity with averaging 
procedure with a classical probability distribution, the function P{a) is not a true proba- 
bility distribution. Because of the overcompleteness of the coherent states, the P function 
is not a unique, well-behaved positive function for all density operators. 

Cahill [63] studied the P representation for density operators which represent pure 
states and found a narrow class of states for which the P representation exists. These states 
can be generated from a particular coherent state |q; > by the application of a finite number 
of creation operators. For example, for the ground state of the harmonic oscillator it is easy 
to show that Xa^ (^) = 1 for all A. In that case the P function becomes P{a) = 5'^ (a). The 
delta function and its derivatives are examples of a class of generalized functions known 
as tempered distributions [61]. Also Cahill [64] introduced a representation of the density 
operator of the electromagnetic field that is suitable for all density operators and that 
reduces to the coherent state P representation when the latter exists. The representation 
has no singularities. 



D{A, 1) = exp(Aa+) exp(-A*a). 



Xiv(A) = IV[pexp(Aa+)exp(-A*a)] = x(A, 1), 




< a 
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Sudarshan [57] offered a singular formula for the P representation in terms of an 
infinite series of derivatives of the delta function. From the mathematical point of view, 
such a series is usually not considered to be a distribution function [61,62]. 

For s = — 1 we have 

^(A, -1) = exp(-A*a) exp(Aa+). 

The s ordering corresponds to antinormal ordering. Because the Q function is the Fourier 

transform of the characteristic function 

Xa(A) = IV[pexp(-A*a) exp(Aa+)] = x(A, -1), 

it follows from Eq. (3.48) that the distribution ^{a, —1) is the Q function. It is given by 
the diagonal matrix elements of the density operator in terms of coherent states: 

Q{a) — — < a\p\a > . 

TT 

Though for all density operators the Q function is bounded, non-negative and infinitely 
differentiable, it has the disadvantage that not every positive Q function corresponds to 
a positive semidefinite Hermitian density operator. Evaluating moments is only simple in 
the Q representation for antinormally ordered operator products. 

For s = 0, the distribution $(A,0) becomes the Wigner function W. The latter 
function is defined as the Fourier transform of the characteristic function 

Xs(A) = IV[pexp(Aa+ - A*a)] = x(A, 0). 

Because this characteristic function is identical to xl^^O)? we conclude that 0) is 
the Wigner function W{a). Therefore, the Wigner function can be used to calculate 
expectation values of symmetrically ordered operators: 



The symmetrically ordered operators are the arithmetic average of {m+n)\/m\n\ differently 
ordered products of m factors of a""" and n factors of a. An example for m = n = 1 is given 
inEq. (3.50). 
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The Wigner function is a nonsingular, uniformly continuous function of ex. for all 
density operators and may in general assume negative values. It is related to the density 
operator as follows: 

W{a) J d^ATr[exp(A(a+ - a*) - A*(a - a))p]. 

Also it can be obtained from the P representation: 

W{a) = ^J P{P)exp{-2\a- P\^)d^p. 

By using the standard transformations [65-67], the master equation (3.7) can be trans- 
formed into a differential equation for a corresponding c-number distribution. Here we 
apply these methods to derive Fokker-Planck equations for the before mentioned distribu- 
tions: Glauber P, Q and Wigner W distributions. Using the relations 

= [(^ - 1) Y + ^^WK s) = D(A, + 1)^ + a+], 

= [(. + 1)| - a]D{A, s) = D{A, s)[{s - 1)| - a], 

one can derive the following rules for transforming the master equation (3.7) into Fokker- 
Planck equations in the Glauber P{s = 1), Q{s = —1) and Wigner W{s — 0) representa- 
tions: 

s — 1 ^ s ~\~ 1 d 

^ ^ 2 da*' ^ ^ 2 da' 

. s+1 d , s-1 d 

pa<r^(a —- — )$, pa^ (a t^^)^- 

^ ^ 2 da*' ^ ^ 2 da' 

Applying these operator correspondences (repeatedly, if necessary), we find the following 

Fokker-Planck equations for the distributions $(«, s) [19]: 

d^(a s) d d 1 d'^ d'^ d'^ 

(3.51) 

Here, ^{a,s) is P{s = l),Q(s = —1) and W{s = 0.) While the drift coefficients are the 
same for the three distributions, the diffusion coefficients are different: 

da = -(A -I- iu})a + iia*, D^a = Di-\- s/i, Daa* = D2 - sX. 
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The Fokker-Planck equation (3.51) can also be written in terms of real coordinates xi and 
X2 defined by 



TflUJ 1 

a = xi + ix2 = \l -r^ < Q > +i = < P >, a* = xi — 1x2 (3.52). 



as follows: 



(3.53) 

with the new drift and diffusion coefficients given by 



di = —{\ — IJ,)xi + UJX2, d2 = —OJXi - (A + IJ,)X2, 

Dif = Uu;D,, - J(A - /x), D!^i = ^ ^ - J (A + ;.), = U 



We note that the diffusion matrix 

r,(s) r)(«) 
-^12 -^22 

for the P distribution (s = 1) needs not to be positive definite. 

Since the drift coefficients are linear in the variables xi and X2{i = 1, 2): 

d --\^A-x- A - ^"^^ 



1 ^^j' 



A= [ ) (3.54) 



with 

u A + /i 

and the diffusion coefficients are constant with respect to xi and X2, Eq. (3.53) describes 
an Ornstein-Uhlenbeck process [68,69]. 

The solution of the Fokker-Planck equation (3.53) can immediately be written down 
provided that the diffusion matrix D is positive definite. However, the diffusion matrix in 
the Glauber P representation is not, in general, positive definite. For example, if 

DriD22 - (0(2)' < 0, 



the P distribution does not exist as a well-behaved function. In this situation, the so-called 
generalized P distributions can be taken which are well-behaved, normal ordering functions 
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[51]. The generalized P distributions are nondiagonal expansions of the density operator 
in terms of coherent states projection operators. The Q and W distributions always exist; 
they are Gaussian functions if they are initially of Gaussian type. 

From Eq. (3.53) one can directly derive the equations of motion for the expectation 
values of the variables xi and X2{i = 1,2): 

, 2 

=-^A,,<a;,>. (3.55) 

By using Eqs. (3.5), (3.52) and (3.55) we obtain the equations of motion for the expectation 
values (Tq{t), o-p{t) of the coordinate and momentum of the harmonic oscillator which are 
identical with those derived in the preceding two subsections by using the Heisenberg 
representation and the method of characteristic function, respectively (see Eqs. (3.13)). 
The variances of the variables xi and X2 are defined by the expectation values 

— Xj^x j X -i X j '^'i 3 — 1 5 2 . 

They are connected with the variances and covariance of the coordinate q and momentum 
P by 

aqq = {2h/mu;)aii, app = 2hmoja22, CTpq =< ^{pq + qp) > - <P >< q >= 2hai2 

and can be calculated with the help of the variances of the quasiprobability distributions 
(*,i = l,2): 

a\^^ = J XiXj^{xi,X2,s)dxidx2 — J Xi^{xi,X2,s)dxidx2 J Xj^{xl,X2,s)dxldx2■ 
Th.e following relations exist between the various variances: 



(^ii - (^ii + ^ - (^ii - ^ - (^ii 1^ - 1>2, (712 - (^12 - (^12 - (^12- 



The variances cr^J-* fulfill the following equations of motion: 

da^'^ 



dt 

1=1 
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They can be written explicitly in the form: 



— —2A2ia[^2 ~ 2yl22C22'' + -^22^ 



dt 

— -1^11 + ^22j'7i2 -^2lCrii - ^12<722 +-^12' 

where the matrix elements are defined in Eq. (3.54). These relations are sufficient to 
prove that the equations of motion for the variances an and (722 and the covariance ai2 are 
the same for all representations as expected. The corresponding equations of motion of the 
variances and covariance of the coordinate and momentum coincide with those obtained 
in the preceding two subsections by using the Heisenberg representation and the method 
of characteristic function, respectively (see Eqs. (3.14)). 

In order that the system approaches a steady state, the condition X> v must be met. 
Thus the steady state solutions are 

^{xi,X2,s) = - — , exp[-- ^ {a-^)ij{oo)xiXj], (3.57) 

27rVdet(or(oo)) ^ . j^^^^ 

where the stationary covariance matrix 



12 (OO) 
22 



^^^^(oo) 



can be determined from the algebraic equation (see Eq. (3.56)): 



1=1 



With the matrix elements Aij given by (3.54), we obtain 

(s). , _ (2A(A + /x) + a>^)Dff + u^Pjl^ + 2a;(A + /x)L>g 
•^11 1°°/ — 



4A(A2 + a;2 - //2) 

'22 (oo) = 



(s) _ oj'^d[1^ + (2A(A - //) + uj'')Dil^ - 2uj{X - ii)d[1^ 



4A(A2 + a;2 - /j,^) 



-a;(A + fi)D[f + u{X - /x)Dg + 2(A2 - ^2)^;^ 
■^2^""^" 4A(A2 + a;2-/x2) 
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The explicit matrix elements cr^j'^ for the three representations P, Q and W can be obtained 
by inserting the corresponding diffusion coefficients. The distribution functions (3.57) can 
be used to calculate the expectation values of the coordinate and momentum and the 
variances by direct integration. The following relations are noticed [48]: 

3.4 The master equation in the Weyl-Wigner-Moyal representation 

The Weyl-Wigner-Moyal representation is a remarkable phase-space representation of 
the quantum mechanics. Roughly speaking, a phase-space representation of the quantum 
mechanics is a mapping from the Hilbert space operators to the functions on the classical 
phase-space which is such that if A is mapped onto /a (a;, y) and p is mapped onto fp{x, y), 
then 

/•oo POO 
/ fp{x,y)fA{x,y)dxdy. (3.58) 
-oo J — oo 

In reality, it is not exactly so, because the Weyl mapping is a mapping from the functions 
on the phase-space to the Hilbert space operators. This mapping W was defined by Weyl 
[70] in the following way: 

-| POO PCO POO POD 

^(/) = 7^r^/ / W / e~^'/^^(-^-y^^f{x,y)dxdy)W{C,r])dCdr]. 
From this it follows very formally that for any p G V{7i), 

1 poo POO POO poo 

T^(P^(/)) = T^TW / / f(^^y)( / e-^'/''^^-^-y^^Tr{pW{^,v))d^dv)dxdy 

(3.59) 

and that (3.59) can be put in standard form (3.58) if the following function on the phase 
space is associated to any p e ViTi), 

1 pOO pOO 

/p(^,2/) = ^^y_ j_ e-(^/^)(-''-^^)TV(pW^(e,r7))derfr7. (3.60) 

The mapping p ^ fp defined by (3.60) is exactly the Wigner mapping [54] which is the 
dual of the Weyl mapping / — > W{f) (hence it can be denoted by W) and fp = W{p) 
is the Wigner function corresponding to the quantum state p e 'D{7i). But the quantum 
nature of the expectation value is not lost because fp{x, y) is not a probability distribution 
on the phase space, taking positive and negative values. 
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In the following, the phase space representation of the master equation (3.3) is ob- 
tained [17] by using the Wigner mapping (3.60). Denoting 



fix,y,t) = fp(t){x,y) = U,(p){x,y), 
it follows from the definition (3.60) that 

-1 fOO /•OO 

/(^,Z/,t) = ^^y_ J_ e-^'/^^^^^-y^^Tripit)Wi^,rjmdv (3-61) 



Then 

dfix,y,t) ^ '■^ '■^ 
dt 

and by duality 



-1 nOO /"OO 



^ (2^ /I /I e-(^/^)(-''-^^)lV(p(t)L(W^(e, v)mdv. (3.62) 

But from (3.10) in the case t = 0,A = W{^, rj) and by using the well-known identities for 
the Fourier transformation, Eq. (3.62) is transformed into the following evolution equation 
for the Wigner function [17]: 

dfix,y,t) y df{x,y,t) ^ df{x,y,t) ^ ^ / ^/ +^^ 

+(^ + ^'^)) + + + ^D,, g^g^ ■ (3.63) 

This equation looks classical and in fact is exactly an equation of the Fokker-Planck type. 
But not every function f{x,y,0) on the phase-space is the Wigner transform of a density 
operator. Hence, the quantum mechanics appears now in the restrictions imposed by 
this last condition on the initial condition f{x,y,0) for Eq. (3.63). Unfortunately, these 
restrictions are not known explicitly. 

Because the most frequently used choice for /(a;,y,0) is a Gaussian function and 
because Eq. (3.63) preserves this Gaussian type, i.e., f{x, y, t) is also a Gaussian function, 
the differences between the quantum mechanics and classical mechanics are completely lost 
in this representation of the master equation. This is an explanation for the frequently 
occurring ambiguities on this subject in the literature. 
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Now from (3.61) by duality it follows that 
and from the results of Subsection (3.1) it follows that 

-| noo poo 



where 

m = am + mv, vit) = im + mv, 

If the initial state p(0) is a pure state corresponding to a coherent wave function 



p(O)0 = ('i/'a,(0),<7p(0), </>)'0a,(O),<7p(O) 

centered m. x — aq{Q)^y — crp(O), i.e., if 
where 

i;o{x) = (27ra,,(0))-i/^exp(-xV4a,,(0)), 

the Wigner function can be analytically calculated because the integrand is an exponential 
function with the exponent having a quadratic form in ^ and rj. 

With the help of the Wigner function the coordinate and momentum probability 
distribution are defined, respectively, by 

/oo 
f{x,y,t)dy 
-oo 

/oo 
f{x,y,t)dx. 
-oo 

^/ X 1 / (x - aa(t))^ . 

Pix,t) = exp — tH^) 

^ ' v/27ra,,(t) 2aqq{t) > 
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and 



It follows 



and 

^/27^app{t) 2(Tpp{t) 
where, with the notations defined in Subsection (3.1) 

agg{t) = S\t)agg{0) + /^^ (t ) (7pp (O ) + 2B{t), 

(7pp{t) = a^{t)(7pp{0) + ^^{t)agg{0) + 2A{t) 

and o"pp(0) — ti^ /4aqq{0). With these notations the Wigner function takes the following 
form [17]: 

= / ^ exp{- ] — 

2nJapp{t)aqq{t) - alq{t) 2{app{t)agg{t) - ^^.(t)) 



x{app{t){x - aq{t)f + aqq{y - ap{t)f - 2apq{t){x - aq{t)){y - cTp(t)))}, 

where 

f^pg(^) = + cx{t)(3{t)app{Q) + 2C{t). 

Evidently 

/oo /"OO 
/ f{x,y,t)dxdy = 1. 
-oo J —oo 

3.5 Density matrix of the damped harmonic oscillator 

In this Subsection we explore the general results that follow from the master equation 
of the one-dimensional damped harmonic oscillator. Namely, we discuss the physically 
relevant solutions of the master equation, by using the method of the generating function. 
In particular, we provide extended solutions (including both diagonal and off-diagonal 
matrix elements) for different initial conditions [20]. 

The method used in the following follows closely the procedure of Jang [21]. Let us 
first rewrite the master equation (3.7) for the density matrix by means of the number 
representation. SpecificalUy, we take the matrix elements of each term between different 
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number states denoted by \n >, and using a+jn >= ^/n^\^\n+l > and a\n >= ^/n\n—l >, 
we get 

= -iujim - n)pmn + >^Pmn - (m + 77, + l)D2pmn 

--^Jm + l)nDlpm+i,n-i - •\/m(nTI)-DiPm-i,n+i 
+ ^V(n + l)(n + 2)(L'i - /i)p^,n+2 + {m + l){m + 2){Dl - p)pm+2,n 

+ - + P)pm-2,n + ^\/(n- l)n(i:>i + IJi)pm,n-2 

+ ^J {m + l){n + 1){D2 + A)p^+i,„+i + ^/nm{D2 - X)pm-i,n-i- (3.64) 

Here, we have used the abbreviated notation =< m\p{t)\n > . This master equation is 
comphcated in form and in the indices involved. It comprises not only the density matrix 
in symmetrical forms, such as Pm±i,n±i, but also matrix elements in asymmetrical forms 
like Pm±2,nj Pm,n±2 and Pm:fi,n±i- In Order to solve Eq. (3.64) we use the method of a 
generating function which allows us to eliminate the variety of indices m and n implicated 
in the equation. When we define the double-fold generating function by 

G{x, y,t)^Yl ^i=^"^l/>mn(t), (3.65) 
— Vmlnl 

the density matrix can be evaluated from the inverse relation of Eq. (3.65): 

Pmnit) = -yl={-^r{^rGi^,y,t)\.=y=0, (3.66) 

Vmlnl ox dy 

provided that the generating function is calculated beforehand. When we multiply both 



sides of Eq. (3.64) by x"^y^/ \J mini and sum over the result, we get the following linear 
partial differential equation of second order for G{x, y, t): 

-Gix, y, t) = {[-{lu + D2)x - Dly] — + [-Dix + {lu - D2)y] — + {D2 + X)- 



dt dx dy dxdy 

+ \^^^*^-^')^ + ^^^-^')^W\{Dl+p)x^ + \{Dl+p)y^ 

(3.67) 

The form of the solution of the Fokker-Planck equations treated in [17] and [21] suggest 
us to look for the solution of Eq. (3.67) of the form 

G{x, y,t) = ^ exp{a;y - [B{x - Cf + D{y - Ef + F{x - C){y - E)]/H}, (3.68) 
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where A, B,C, D, E, F and H are unknown functions of time which have to be determined. 
When we first substitute the expression (3.68) for G{x,y,t) into Eq. (3.67) and equate 
the coefficients of equal powers of x, y and xy on both sides of the equation, we get the 
following differential equations for the functions A,B,D,F and H: 

-J^- -(Dl -/.)§- (I)i -/.)§- (1)2 + A)^ + 2A, (3.69) 

|(|) = 2(A -iu;)^-„^- ^-{D, - - 2{D, + A)^ - 2(iP* - (3.70) 

-(-) = 2(A + iu;)--i.-- -{Dl 2{D, + A)-^ - 2{D^ -/.) — , (3.71) 

rl TP TP TP'^ fi ]~) J~) F{ HP J~) P 

J,iji) = 2A--(fl,+A)^-2M^ + ^)-4(Z>.+A)^-2(DI-,)^-2(D,-,)^. 

(3.72) 

In addition to these equations, we get for the functions C and E 

2B— + F— = (-2(A - iu;)B + iiF)C + (2//S - (A + iu;)F)E, (3.73) 

(JjL (JjL 

+ F^ = (-2(A + iu;)D + iiF)E + (2ij,D - (A - iu;)F)C. (3.74) 

ti/i' (J/V 

The equations (3.73) and (3.74) can be reformulated in order to eliminate the functions 
B, D and F, provided BD - F^/4 ^ 0. We obtain 

^ = -(A-iw)C + //E, (3.75) 
dE 

= -(X + iu;)E + iiC. (3.76) 

(JjL 

The functions A,B,D,F and H are connected by the auxiliary condition that Trp is 
independent of time. The trace of p can be evaluated by summing the diagonal matrix 
elements pnn given in Eq.(3.66) or directly by using the integral expression 

CX) If 

Trp = Pnn = TTT-To / exp(-A;iA;2) exp(iA;ia; + ik2y)G{x, y, t)dkidk2dxdy. 
t^o (27r) J 

We obtain with the generating function (3.68) 

442 p2 

'^P={^{^-BD))-^/\ (3.77) 
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This quantity is time-independent which can be verified by constructing an equation sat- 
isfied by the quantity {F'^/A - BD)/H'^. Combining Eqs. (3.70)-(3.72) we get 

We see immediately that the first factor on the right-hand side of this equation is identical 
with the right-hand side of Eq. (3.69). Accordingly, we find 

Since the scaling function H is arbitrary, we simplify the following equations by the choice 

— -BD = -H. (3.78) 

Setting Trp = 1, we obtain from Eqs. (3.77) and (3.78) the normalization constant A'^ = 
—H/4. As a consequence, we can simplify Eqs. (3.70)-(3.72) by eliminating the function 
H from these equations. The resulting three equations are 

-2(A + iio)B -nF + 2{Di - //), 

■2{X-iu)D - i^F + 2{Dl- 1^), (3.79) 

2^{B + D)- 2XF - A{D2 + A). 

These equations imply that the function D is complex conjugate to B, provided that the 
function F is real. 

In order to integrate the equations for the time-dependent functions B,C,D,E and 
F we start with Eqs. (3.75) and (3.76). These equations imply that the function E is 
complex conjugate to the function C. By solving the coupled equations we find: 

C{t) = E*{t) = u{t)C{0) - v{t)C*{0), (3.80) 

where u{t) and v{t) are given by (3.45) and (3.46) for the two considered cases: overdamped 
and underdamped, respectively. For integrating the system (3.79) we proceed in the same 
way as for integrating the system (3.14). With the assumption that F is real and 



It 

dD 

It 
dF 
'dt 



D{t) = B*{t) = R{t) + il{t), 
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we obtain explicitly: 



272^ 2 ^ 27 



F{t) = -e-''^'^{2^iI + uF) 

7 



7^ 27 
where we used the notations: 



At± = A ± 7, 7 = ^J^F^^, R = i?(0) - i?(oo), / = /(O) - /(oo), F = F(0) - F(oo). 

We can also obtain the connection between the asymptotic values of B{t), D{t), F{t) and 
the coefficients Di, D2, A* ^^id A: 

A(ReDi - /x) + a;Imi:>i + //(Da + A) 



R(oo) = ReD(oo) = 

/(oo) = ImL>(oo) = 
F(oo) = -2 



A2 - 72 

a;A(Re/:»i - //) + (//^ - X'^)ImDi + 0Jn{D2 + A) 
A(A2-72) ^ 

/i[A(Re/:>i - //) + colmDi] + [X^ + uj^){D2 + A) 



A(A2-72) 

When all explicit expressions for A, B, C, D, E, F and H are introduced into Eq. (3.68), we 
obtain an analytical form of the generating function G{x, y, t) which allows us to evaluate 
the density matrix. 

If the constants involved in the generating function satisfy the relations 

C(0) = 0, R(0) = R(oo), 1(0) = /(oo), F(0) = F(oo), 

we obtain the stationary solution 

C{t) = E{t) = 0, R{t) = R{0), I{t) = /(O), F{t) = F(0) 

and 

D(t) = B*{t) = R{0) + z/(0), H(t) = -AA^(t) = R^(0) + /^(O) - F^(0)/4. 
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Then the stationary solution of Eq. (3.67) is 

1 r,. 1 



Gix, y,t) = j exp{(l - ^)xy - {Bx^ + B*y^)}. (3.81) 
In addition, for a thermal bath [17] with 

rruvDgq = ^12- Dpq = 0, = 0, 
the stationary generating function is simply given by 

The same generating function can be found for large times, if the asymptotic state is a 
Gibbs state with = 0. In this case we obtain with Eq. (3.26) and = 

Do ~ A coth — — 
2kT 

and 



G{x, 2/) = (1 - exp(-^)) exp(exp(--^)a;y). 

The density matrix can be calculated with Eq. (3.66) and yields the Bose-Einstein distri- 
bution 

Pmn(oo) = (1 - exp(- — )) exp( rj^)5mn- 

A formula for the density matrix can be written down by applying the relation (3.66) to 
the generating function (3.68). We get 

Pmn{t) = exp[-(SC2 + DE' + FCE)/H] 

^ _ F\n:^f B \m r D\n-:,r 2BC , \m-2ni-n3 / , FC_\n-2n2-n3 

^ „ ni!n2!n3!(m — 2ni — n3)!(n — 2n2 — ns)! 

ni,n2,n3=0 

(3.82) 

In the case that the functions C{t) and E{t) vanish, the generating function has the form 
of Eq. (3.81). Then the elements of the density matrix with an odd sum m + n are zero: 
Pmn = for m + n = 2/c + 1 with k = 0, 1, 2, ... The lowest non- vanishing elements are 
given with pmn = Pnm as 

POO -J, P20--^, Pn-]4(1-^)' 
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2BB* 1 , F 2 F ^/6B ^/6B^ 

P22 - + , P31 - -(1 - Jj)^, P40 - 

It is also possible to choose the constants in such a way that the functions B and D vanish 
at time t = and F(0) = H{0). Then the density matrix (3.82) becomes at t = 0{E = C*): 

Pmn(O) = ^i=(C*(0))-(C(0))"exp(-|C(0)|2). (3.83) 

This is the initial Glauber packet. The diagonal matrix elements of Eq. (3.83) represent a 
Poisson distribution used also in the study of multi-phonon excitations in nuclear physics. 
In the particular case when we assume 

Di=n = 0, ^2 = A, B{0) = D(0) = 0, F(0) = H{0) = -4, 

the differential equations (3.79) yield B{t) = D(t) = and F{t) = H{t) = -4. Then the 
density matrix subject to the initial Glauber packet is (see also [21]) 

pmn{t) = -^{c*{t)r{c{t)reM-\c{t)n 

where C(t) is given by Eq. (3.80). 

4. Applications to nuclear equilibration processes 
4.1 Charge equilibration in deep inelastic reactions 

Over many years experimental data have been measured in the field of deep inelastic 
heavy ion collisions [71,72]. The characteristic feature of these collisions is the binary 
character of the system, i.e. the final fragments have nearly the same masses as the initial 
nuclei. Deep inelastic collisions are mainly described by the dynamics of selected collective 
degrees of freedom. These are the relative motion of the nuclei, mass and charge exchange, 
the neck degree of freedom and surface vibrations of the fragments [73] . Another important 
feature of these reactions is the dissipation of energy and angular momentum out of the 
collective degrees of freedom into the intrinsic or single-particle degrees of freedom. 

In most of the approaches the collective degrees of freedom are chosen in an a 'priori 
way, guided by measured macroscopic variables which do not cover the complete space 
associated with all of them. Also, intrinsic degrees of freedom play a role. Thus, by taking 
into account the coupling between the actually treated and intrinsic degrees of freedom 
and other collective degrees of freedom, irreversible processes occur in the dynamics of the 
actually treated collective variables. 
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One of the most characteristic processes in deep inelastic reactions is the fast redis- 
tribution of neutrons and protons among the coUiding nuclei in the early stage of the 
reaction called charge equilibration, neutron excess mode or isospin relaxation. The con- 
ditional variance of the charge distribution at fixed mass asymmetry has been largely used 
as a sensitive quantity to investigate the nature of the charge equilibration process. 

There are two opposite approaches which have been developed in order to describe 
the charge equilibration mode, namely the quantum mechanical collective treatment and 
transport theories. The collective treatment stresses the coherent properties of the pro- 
cess and considers this mode as connected to the isovector dipole giant resonance of the 
composite system. An experimental indication for a quantum process is that the charge 
variance reaches a saturation value at large values of the total kinetic energy loss. The 
other opposite view relies on stochastic exchange processes between the colliding heavy 
nuclei. 

In the following the charge equilibration mode is treated as a one-dimensional collective 
mode in the charge asymmetry degree of freedom. The damping of this collective mode is 
due to the coupling to the remaining collective and intrinsic degrees of freedom. 

One method of introducing dissipation in a quantum mechanical description of deep 
inelastic collisions is to assume that the energy dissipation is similar to the loss of energy 
of a harmonic oscillator coupled with a large number of other harmonic oscillators. This 
mechanism can be simulated by a friction term of the Kostin type in the Schrodinger 
equation. According to this line, the role of collective motion and quantum fluctuations 
for charge and mass equilibration in deep inelastic reactions has been investigated [74,75]. 

Another method of treating dissipation in quantum mechanics, especially for the de- 
scription of coupled collective modes, is the axiomatic method of Lindblad. In this method, 
as we said before, a simple dynamics for the subsystem of the explicitly treated collective 
degrees of freedom is chosen, namely a semigroup of transformations which introduces a 
preferred direction in time and, therefore, can describe an irreversible process. 

In the present Subsection, Lindblad's theory for open quantum systems is applied 
to the problem of the damping of a single collective degree of freedom in deep inelastic 
collisions. We show the description of the charge equilibration mode by a damped harmonic 
oscillator in the framework of the Lindblad theory [17,18,22]. 
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For a comparison of the theoretical results with experimental data we consider the 
charge equilibration in deep inelastic collisions of heavy ions. We can take over the solutions 
for the centroids and variances of a damped oscillator given by (3.15), (3.16), (3.19), (3.27), 
(3.28) for a comparison in order to illustrate our results. In [22] the charge equilibration 
of the systems ^^Fe +^^^ Bi,^^ Fe +^^^ U and ^^Mo Mo was analysed within a collective 
description of the charge asymmetry degree of freedom. The charge asymmetry coordinate 
is defined as 



where Z\ and Z-i are the charge numbers of the fragments. The Hamiltonian used for the 
description of this degree of freedom is a harmonic oscillator Hamiltonian of the form 



The stiffness parameter K and mass B were calculated in [22] on the basis of the liquid 
drop model for the potential energy and by use of a semiempirical formula for the frequency 
u) = {K/By/'^. For the system ^^Fe +2°^ Bi the following values were obtained [22]: B = 
1.6127 X 10~^° MeV s^ and K = 1.6363 x 10^ MeV, corresponding to an oscillator energy 
fioj = 6.63 MeV. Figs. 1-3 show the charge centroids and charge asymmetry variances as 
functions of time for the reaction of ^^Fe on ^°^Bi at Eiab = 465 MeV. The full dots are 
the experimental data of Breuer et al. [76] . Since the experimental data were obtained as 
functions of the total energy loss. Pop et al. [22] have related the data with the reaction 
times in order to compare them with calculations. As a model they used a parametrization 
of the deflection function in terms of relaxation times for the dissipation of the radial kinetic 
energy and relative angular momentum and for the development of deformations [77-79]. 
These times were fixed by reproducing the ridge of the double-differential cross-section 
(Pa/dQdE plotted as a contour diagram in the total kinetic energy against scattering 
angle (Wilczynski plot). For short reaction times the kinetic energy loss is proportional 
to the reaction time (~ 30 Mev/10~^^s). Our comparison does not depend much on the 
specific assumptions of the model used to derive the reaction times from the measured 
total energy losses. 

The "experimental" values of the variance in Fig. 3 are constructed from the experi- 



Zi — Z2 

ZI + Z2 
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mental values of cr^^ by using the equation (see Eqs. (3.14)) [22] 



2(A - ij,)a^^ + -a^p^ + 2D^^ 



dt 



and applying the corresponding fitted values for the parameters A, and D^^. The curves 
are calculated with the parameters ^o,A and fi fitted to the experimental data in Fig. 1 
and in addition with u>i = —2mujDqq/Ti,u}2 = —2Dpp/m?ko and uj^ = 2Dqp/h fitted to the 
experimental variances in Figs. 2 and 3 (here q = $,tP — p^Tiri — B). The set of parameters 
U3i is restricted by the condition (3.4) which now reads as 



The full curves in Figs. 1-3 represent the overdamped solution (a; = 1.0073 x lO^^s"^ < 
H = 2.33 X lO^^s"^) for the centroids and variances. This solution is the only one for 
which a good description of all the experimental data could be obtained. In the case of the 
underdamped solution no set of parameters could be found fitting simultaneously the data 
for the centroids and variances. Similar results were also obtained for the other considered 



The parameters of the theory are accepted as free quantities, fitted to the experimental 
data of the charge distributions. The overdamped solution with imposed fundamental 
constraints for the damped quantum oscillator succeeded to describe the experimental 
data within quite restricted ranges of the model parameters. Both the underdamped and 
the equilibrium solutions can not describe the data for the charge centroids and the isobaric 
variances. 

The given treatment has the advantage that analytical solutions can be obtained 
for the damping of a harmonic oscillator. This is in contrast to the method when one 
introduces a friction term of the Kostin type into the Schrodinger equation and solves 
the resulting equation numerically. However, future investigations have to be carried out 
to interpret the fitted parameters of the theory from a microscopical point of view. For 
example, the parameters can be related to those obtained from a microscopical derivation 
of Fokker-Planck equations. 

4.2 Charge and neutron equilibration and the open quantum system of two 
harmonic oscillators 



UJ1UJ2 — C1J3 > A^. 



systems 



56pg ^238 u and ^^Mo Mo in [22]. 
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In the present Subsection we extend the work of Gupta et al. [80] on the dynamics 
of the charge equihbration process in deep inelastic coUisions and treat the damping of 
the proton and neutron asymmetry degrees of freedom with the method of Lindblad [23] . 
The charge and mass distribution in di-nuclear systems can be described with continuous 
coordinates of the charge and neutron asymmetries defined by 

Zi - Z2 iVi - N2 

Here, Zi{Ni) and Z2(A^2) are the total charge numbers (neutron numbers) on the left-hand 
and right-hand side of a plane through the neck of the di-nuclear system [73]. Without 
damping, the charge and neutron asymmetry degrees of freedom are described by a wave- 
function 'ilj{r]z,'r]N,t), which is the solution of the following time-dependent Schrodinger 
equation: 

H{'qz,m)i'{^z,r]N-,t) = ihdil;{'r]z,r]N,t)/dt, 
where the Hamiltonian of the model is given as {\kzN\ < {kzkN)^^^) '■ 

Hivz.VN) = - + l^^^z + IkWN - kzNVzVN. (4.1) 

The Hamiltonian has the simple structure of two coupled oscillators in the coordinates r]z 
and ?7iv in order to keep the time development of the wavefunction analytically solvable. 
The string constants of the potential can be calculated with the liquid drop model for the 
sticking configuration of the nuclei, i.e. for the relative distance R = Ri + R2, where Ri 
and i?2 are the radii of the two colliding nuclei (for more details see [80]). 

Hofmann et al. [39] have considered the coupling of the mass asymmetry coordinate 
to the charge asymmetry coordinate on the basis of a density operator formalism using a 
quantum master equation in a perturbative treatment. The coupling of the neutron and 
charge asymmetry coordinates has also been studied in the framework of a two-dimensional 
Fokker-Planck equation by Gross and Hartmann [81], Schroder et al. [82], Birkelund et al. 
[83] and Merchant and Norenberg [84]. 

With the Hamiltonian (4.1) as an example, we study the damping of two coupled 
oscillators in the framework of the Lindblad theory. In order to have a formalism which is 
generally applicable, we give the following formulae in terms of the two general coordinates 
qi and q2 instead of rjz and 77 jv- We present the equation of motion of the open quantum 
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system of two oscillators in the Heisenberg picture. With this equation we derive the time 
dependence of the expectation values of the coordinates and momenta and their variances. 
The connection with the Wigner function and Weyl operator is also discussed. Finally, we 
demonstrate the time dependence of the various quantities for a simplified version of the 
model, where the decay costants can be calculated analytically. 

If is the dynamical semigroup describing the time evolution of the open quantum 
system in the Heisenberg picture, then the master equation is given for an operator A as 
follows [14,16,17]: 

^^-^^ = ^^^'^^^^ = + ^Y.(Ki^M),Vj\ + [V+MA)]V,). (4.2) 
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The operators H, V, , Vj^{j = 1, 2, 3, 4) are taken to be functions of the basic observ- 
ables of the two quantum oscillators. The coordinates qk and the momenta pk obey the 
usual commutation relations {k,l = 1,2): 

[Qk,pi] = ifi^ki, [qk, qi] = \pk,pi] = 0. 

In order to obtain an analytically solvable model, H is taken to be a polynomial of second 
degree in these basic observables and Vj,V^ are taken to be polynomials of only first 
degree. Then in the linear space spanned by qkiPk{k = 1,2), there exist four linearly 
independent operators 1^=1,2,3,4: 

2 2 

^ X] '^jkPk + ^ bjkqk, 
k=i fc=i 

where Oj-fc, 6^7^ G C with j = 1, 2, 3, 4. It yields 

2 2 

fe=i fe=i 

where a^^, 6*^ are the complex conjugates of ajk, bjk- 

The Hamiltonian H is chosen in the form of two coupled oscillators 

2 2 1 2 

TT \ — ^/ -'- 2 '^k^h 2\ 7 -'- \ — ^ / \ 

^ = 2^(^^Pfe + ^— ?fc) + ^i2PiP2+ 2 ^ikukAPk^qk2 + qk^Pki) + J^i2qiq2- 

k=l ^ fei,fc2 = l 
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Here we use the following abbreviations {k = 1,2) : 



D 



-Re(afeai), 



Dr 



h 



Re{hlhi), D, 



h 



qkQi ~ ^qiQk ~ 2 ^ '•'^ ^PkPi ~ ^piPk ~ 2 ^ '■'^ ^ikPi 

ai2 = -Q!2i = -Im(aia2), P12 = -P21 = -Im(bib2), Xki = -Im(afcb/). (4.3) 

The scalar products are formed with the vectors a^jb^ and their complex conjugates 
a^, b^. The vectors have the components 

afc = (aifc, a2fc, osfc, a^k), bfc = (bik, ^2fe, hk, b4k). 

Now, as a consequence of the definitions (4.3) of the phenomenological constants which 
appear in L{A) and of the positivity of the matrix formed by the four vectors ai, a2, bi, b2, 
it follows that the principal minors of this matrix are positive or zero. This matrix is given 
by: 

/alai a^a2 a^bi a^'ba \ 



D 



D 



a^ai a^as a^bi a^ba 

hlai b^a2 b^bi b];b2 

Vb^ai b^a2 b^bi b^b2/ 

Dqiq2 - ihai2/'2 -Dq^p^ - Z^Aii/2 --Dgip2 - ihXi2/2\ 



92 91 



D 



D 



piqi 



iha2i/2 
ihXii/2 -Dp^q^+inX2i/2 



q2q2 



-Dq,p, - ihX2i/2 -D. 



q2P2 



ihX22/2 



D 



\ -Dp^q^ + inXi2/2 -Dp^q^ + ihX22/2 Dp^p^ - ifi/32i/2 



Dp^p^ - ihPi2/2 



D. 



P2P2 



J 
(4.4) 



For example, we can write one of the conditions obtained from the positivity of (4.4): 

-^9191-^92 92 ~ -^9192 — 4^ '^12- 

This inequality and the corresponding ones derived from (4.4) are constraints imposed on 
the phenomenological constants by the fact that $t is a dynamical semigroup [14,16,17]. 

The time-dependent expectation values of self-adjoint operators A and B can be writ- 
ten with the density operator p, describing the initial state of the quantum system, as 
follows: 

ruAit) = Tr{p^t{A)), aAsit) = ^Tr{p^t{AB + BA)). 

In the following we denote the vector with the four components mg. (t), m^. {t),i — 1,2, 
by ni(t) and the following 4x4 matrix by (j{t): 



( 



(T. 



ait) 



qiqi 
'^9291 



(T, 



9i92 
92 92 



'9lPl 
''92P1 



'91P2 
''92P2 



^Piqi '^pi92 "^pipi 



P291 



'P292 
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'P2P1 



^PlP2 
'^P2P2 



Then via direct calculation of L(qk) and L(pk) we obtain 

dm 



dt 



Ym, (4.5) 



where 

-A21+//2I — A22+A*22 «12 + ^12 1/"T.2 

-miuf P12 - 1^12 -All - A*ii -A21 - At2i 

V -P12 - 1^12 -m2ijJ2 -A12 - Atl2 -A22 - At22 . 



Y = 



(4.6) 



From (4.5) it follows that 

m{t) = M(t)m(0) = exp(tr)m(0), (4.7) 

where iii(O) is given by the initial conditions. The matrix M{t) has to fulfil the condition: 

lim M{t) = 0. (4.8) 

In order that this limit exists, Y must have only eigenvalues with negative real parts. 
By direct calculation of L^q^qi), L{pkPi) and L{q]^pi + Piqk), {k,l = 1, 2), we obtain 

dcr ^^ 

— =Ya + aY^ + 2D, (4.9) 
where D is the matrix of the diffusion coefficients 



D 



and the transposed matrix of Y. The time-dependent solution of (4.9) can be written 

as 

(j(t) =M(t)((j(0)-E)M^(t) + E, (4.10) 

where M{t) is defined in (4.7). The matrix E is time independent and solves the static 
problem (4.9) {da/dt = 0): 

Y'S + Y:Y^ + 2D = 0. (4.11) 
Now we assume that the following limit exists for t — > 00: 

f7(oo) = lim a{t). (4.12) 
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^qiqi 


^qiq2 


^qiPi 


D 




^qiqi 


^q2P\ 


D 


^piqi 


^piq2 


^PiPi 




^P2qi 


^P2q2 


Dp2pi 





In this case it follows from (4.10) and (4.8): 

a(oo) = E. (4.13) 
Inserting (4.13) into (4.10) we obtain the basic equations for our purposes: 

a(t) = M{t){a(0) - a(oo))M^{t) + a{oo), (4.14) 

where 

Ya{oo) + a{oo)Y^ = -2D. (4.15) 

Now we want to discuss the time dependence of the Wigner function. This function is 
defined as: 



-j^ /•OO I'OO I'OO />oo / ^ 



xTT[p^i{W{^i,^2;Vi,V2)Mid^2dvidv2, (4.16) 
where the Weyl operator W is defined by (^1, ^2, ^1, ?72 real) 

VF(^i,^2;?7i,?72) ^ exp[in~^ {ijiqi +77252 - CiPi -6^2)]- 

Using the method developed by Lindblad [14,16,17] for the one-dimensional case, we find 
the following relation for the time development of the Weyl operator: 

^t{W{^u 6; m,V2)) = W{^i{t),Uty,Vi{t),V2{t)) exp^(t). (4.17) 

The real functions ^{t) — {Ci{t), $,2{t); r]i{t) , r]2{t)) and g{t) satisfy the equations of motion: 

di{t)/dt^ JY^J-^i(t) (4.18) 

dg{t)/dt = -n-^at)JDJ-^C{t), (4.19) 

where 

/O -1 \ 
I -1 

1 • 

Vo 1 / 
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Eqs. (4.18) and (4.19) are obtained by inserting the Weyl operator W{^i, ^2',V1jV2) into 
the equation of motion (4.2). The initial conditions for the coordinates ^i{t) , ^2{t) , r]i{t) 
and rj2{t) are determined by ^i(O) = ^i, ^2(0) = ^2, ^i(O) = r]i and r/2(0) = r]2, respectively, 
and g{t) by g{0) = 0. From (4.18) and (4.19) we find that ^{t) is a linear function in the 
coordinates ^i,^2,?7i and 772 and g{t) a quadratic function. 

The Weyl operator can be used to calculate the time-dependent expectation values 
m(t) and a{t) (see (4.7) and (4.10)), since this operator is connected with the coordinates 
and momenta via the derivatives 



dW. 


i 






dW. 


i 


dr]i 






1 




n 



(4.20) 



{PiQj + QjPi), 



d^W _ 1 
ap^'^=°" 2F 

k=0 - -T2^i^3 



driidrjj % 
For example, one obtains by using (4.20) 



Equations of this type can be evaluated with the help of (4.17)-(4.19) and lead to the 
same results for m(t) and a{t) as given before. With the Weyl operator (4.17) we can 
calculate the time development of the Wigner function. For this purpose we use the 
Fourier transform of the Wigner function at t = 0: 

Tr{peKp[ih~'^{r][qi + r/^?2 - C'lPi - C2P2)]} 

/oo pOO i-OO pOO 
/ / / exp[ih~^{xir][ + X2V2 - Vi^i - V^Q] 
-00 J —00 J —00 J —00 

'xfixi,X2,yi,y2,t = 0)dxidx2dyidy2. 



53 



When this relation is inserted into (4.16) after the Weyl operator ^t{W) is expressed by 
(4.17), one can integrate over the coordinates ^i, Vi V2 with the following result for 
the Wigner function: 

2 roG f-oo POO POO 2 

xf(x[,x'2,y[,y!2,t = Qi)dx'^dx'2dy[dy'^, (4.21) 
where x = (a;i, a;2, j/i, ^2) and the matrix Z(t) is given by 

Z{t)= [ M{t')DM^{t')dt'. 
Jo 

This definition can be applied in order to rewrite (4.14): 

a{t) = M{t)a{0)M^{t) + 2Z{t). (4.22) 
In the particular case when we set 

f{xuX2,yuy2,t = Qi) = j^^^^^^^^^^^y^exp[-^(x-m(0))a(0)-\x-m(0))], 
we obtain from (4.21): 

f{xi,X2,yi,y2,t) = [je^(27ra(t))] V2 ^^P[~^(^ ~ " '^(^)]' 

which is the well-known result for Wigner functions [45,69,85]. 

In order to illustrate the developed formalism we present an example of two oscillators, 
which are coupled by a potential of the form as used for the proton and neutron degrees 
of freedom in (4.1), i.e. ki2 = 0, Hij = 0, 1^12 ^ 0. In this case the matrix Y, governing the 
time development of the expectation values m{t) and o"(t), becomes 



Y = 



I -All -A12 1/mi -q;i2 

-A21 -A22 OLx2 1/W2 

-miujl (3i2 - z/12 -All -A21 

V -Pl2 - T^12 -m2C^i -A12 -A22 



(4.23) 



For the calculation of the matrix M{t) we must diagonalize the matrix Y by solving the 
corresponding secular equation, i.e. det(F — zl) = 0, where z is the eigenvalue and / 
is the unit matrix. According to (4.23) one obtains an equation of fourth order for the 
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eigenvalues z, which can be simply solved only for special examples. In the particular case 
with ai2 = 0, Pi2 = 0, Ai2 = and A21 = 0, the secular equation is obtained as 

[(z + Aii)2 + uj] [{z + X22f + ojW = vl^/mim2. 

The roots of this equation have the general structure 

zi = -7+ + za;+, Z2 = -7+ - 10)+, zz = -7- + iu)-, Z4 = -7- - iuj-. 

The constants 7± and u!± can be easily calculated for the case 1/12 = : 

7+ = An, 7_ = A22, UJ+ = uji, a;- = (^2 (4.24) 

or for All = A22 = A,a;i = a;2 = a;(a;^ > i'i2/{'mim2y^^) 

7+ = 7- = A, o;^ = a;^ ± 1^12/ (mim2)^/^. 

Only positive values of 7+ and 7_ fulfil (4.8). Applying the eigenvalues Zi of Y we can 
write the time-dependent matrix M[t) as follows: 

Mmn{t) = J2^rnieM^it)N-\ 

i 

where the matrix N represents the eigenvectors of Y: 

n 

With the relations Mmnit = 0) = 5mn and dMmnit) / dt\t=o = Ymn and using (4.7), (4.14), 
we conclude that the expectation values of the coordinates and momenta decay with the 
exponential factors exp(— 7+^) and exp(— 7_t) and the matrix elements amn with the 
combined factors exp(— 27+t), exp(— 7_t) and exp[— (7+ + j_)t]. 

Since the matrix elements are in general lengthy expressions, we present here 

the matrix M{t) only for the special and simple case that the oscillators are uncoupled. 
With the roots given in (4.24) we obtain 

M{t) = 

/ e-^ii*cosa;it ^— e"-^"* sina;it \ 

— mi^ie"^"* sinwit e~^^^* cosuit 

\ -m2UJ2e~^^^* sinu2t e'-^^^* cosu2t J 

(4.25) 
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This matrix can be used to evaluate a{t) defined by (4.14) or (4.22). For example, we find 
the following expression for cri2 = cr^^g^ with M{t) of (4.25): 

•^9192 (^) = exp[-(Aii + A22)i] ((0-51^2(0) - 0-5,52(00)) cos a;it cos a;2* 

H — (o'g2Pi(0) ~ ^g2Pi(°o)) sina;itcosa;2t H — (<7gip2(0) ~ '^qiP2(°o)) cosa;it sina;2t 

H (crpiP2(0) - o-pj^p2{oo))sinujitsinu2t) + aq^q^{oo). 

fn-l7Jl2UJ\UJ2 

Similar expressions are found for the other matrix elements of cr(t). The matrix 
elements of cr(oo) depend on Y and D and must be evaluated with (4.15) or by the 
relation: 

a{oo) = 2 M{t')DM'^(t')dt'. 
Jo 

As an example we present the value of (Jg^gj 

(00): 

aq^q,{^) = 2{[(An + A22)^ + ((^i + (^2)'] [(An + A22)' + {uJi - uj2)^]}~^ 

x{(Aii + A22)[(Aii + A22)^ +uj^+ ujl]Dq,q^ + [(An + A22)^ + a;? - u;l]Dq^pJmi 
+ [(An + A22)^ +(^2- ^i]-Dqip2/"^2 + 2(An + >^22)Dp^p2/'mim2}. 

Similar expressions are obtained for the other matrix elements of cr(oo). The diffusion 
coefficients Dq^q^^Dq^p^^Dq^p^ and Dp^p^ are in general zero for uncoupled oscillators in- 
teracting with an usual environment. This has the consequence that the expectation values 
'^qiq2-i^qiP2T^q2Pi ^piP2 vanish for t — cxd. It is a very interesting point that the gen- 
eral theory of Lindblad allows couplings via the environment between uncoupled oscillators 
with ki2 = 0, jJLij = 0, 1/12 = 0. According to the definitions of the parameters in terms of 
the vectors a.k and b^, the diffusion coefficients above can be different from zero and can 
simulate an interaction between the "uncoupled" oscillators. In this case a structure of the 
environment is reflected in the motion of the oscillators. 

5. Damping of angular momentum in open quantum systems 
In heavy-ion collisions, especially in deep inelastic reactions, one has experimentally 
studied the loss of angular momentum in the relative motion of the nuclei [72]. Measure- 
ments of the 7-multiplicities yielded information on the spin of the excited nuclei, which 
means information on the loss of angular momentum in the relative motion of the nuclei. 
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If the energy loss in the relative motion is simultaneously measured one can determine a 
time scale of the loss of angular momentum, i.e. one can explore the damping process as 
a function of time. This example motivated us to study the damping of angular momen- 
tum. Another example for the damping of angular momentum is the spin-relaxation in 
the presence of a magnetic field which was theoretically studied in [86]. 

The damping of angular momentum can be described at different levels of approx- 
imations: one may use quantum- mechanical or classical methods and microscopical or 
phenomenological models. In this paper we describe the damping of angular momentum 
with the quantum-mechanical theory of Lindblad [24,25]. We choose the operators, which 
open the system, proportional to the components of the angular-momentum operator, to 
generators of the proper Lorentz group and proportional to a linear combination of mo- 
menta and coordinates. 

Let us assume a finite set of Hermitian operators Ai depending on the basic variables 
of the physical problem and given in the Heisenberg picture. Then their time evolution can 
be determined within the theory of Lindblad by the following Markovian master equations: 

dMt) 



dt 

fe=i 

The matrix elements Mj^ and components Di are time- independent numbers. They are 
calculated by the equations {Ai = Ai{t = 0)) : 

L{Ai) = ^[H,Ai] + ^Y1(^;[A,V,] + [V;,A,]V,) (5.1). 

j 

Here, H is the Hamiltonian, Vj and V^^ are operators which are functions of the basic 
variables. V^^ is the Hermitian adjoint operator to Vj. The expectation values of the 
operators Ai{t) are solutions of the differential equations {i = 1, N) : 

^-^^ = Y^Miu{A,{t)) + Di. 

The following considerations are restricted to three spatial degrees of freedom described by 
the coordinates qi,q2, Qs and linear momenta Pi,P2,P3- We assume that the Hamiltonian 
conserves angular momentum: 

[H,Li]=0, i = 1,2,3. 
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Further we assume for simplicity that the opening of the quantum-mechanical system fulfils 
rotational symmetry. If we restrict this study to three operators ^^=1^2,3, these operators 
have to transform like the components of a vector under spatial rotations. They fulfil the 
commutators 

3 

[5^V+V,-,L,] = 0, i = 1,2,3. 

i=i 

The following three examples for Vj will be considered {j = 1, 2, 3) : 

(1) V, = aLj, 

(2) Vj = aNj = a{spj + qjPs), 

(3) Vj = apj + bqj. 

The operators Nj contain a fourth coordinate s and the corresponding momentum ps = 
—ihd/ds. This coordinate can be thought to describe an intrinsic degree of freedom not 
affected by spatial rotations. The operators Lj and Nj have commutation relations as the 
generators of the proper Lorentz group [87]. The free parameters a, a and b are complex 
numbers. 

(1) For the case Vj proportional to Lj, we obtain from Eq. (5.1) 

L{A) = ^[H,A]+^-^J2^Lj,[A,Lj]]. 

If we set the operator A equal to the square of the angular momentum, ^ = L^, we find 
conservation of the total angular momentum: 

^ = 0. ^ = (8.2) 

dt ' dt ^ ^ 

Here {L'^{t)) denotes the expectation value of L'^{t). The squares of the components of 
get equalized: 

= -\a\'hi3mt) -L% z = 1, 2, 3. (5.3) 
The expectation value of Lf is calculated as 

/ r2\ 

{Ll{t)) = {L'^i)t=oexp{-3h\a\H) + -^{1 - exp{-3h\a\H)), 
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where (L^) is time independent as shown in Eq. (5.2). These results prove that the ansatz 
(1) can be used in apphcations when the square of the angular momentum is conserved, 
but the squares of the individual components are equalized with time. Examples are 
depolarization phenomena. 

We note the fact that in the derivation of Eqs. (5.2) and (5.3) we only used the 
commutators between angular-momentum operators. Therefore, the results in this case 
are also true for systems with many degrees of freedom. 

(2) For the case Vj proportional to Nj we obtain for L{A) according to Eq. (5.1): 

LiA) = ^[H,A] + \^J2^N„[A,N,]]. 

If we choose A = L^, we find the following time derivative of L'^{t) : 

= 2\a\^n{L\t) + N\t)). (5.4) 
Here we used the commutation relations 

3 3 

[Li, Nj] ^ih^ SijiNi, [Ni, Nj] = -ih ^ SijiLi, 
1=1 1=1 

where Siji is the Levi-Civita antisymmetric tensor. Since the right-hand side of Eq. (5.4) 

jj=l 



contains the square N'^ = Yl^j=i '^^ have also to calculate the time derivative of N^{t) 



= ^[H, N\t)] + 2\a\'hiL\t) + N\t)). (5.5) 

In order to simplify the discussion of Eqs. (5.4) and (5.5) let us assume that the commutator 
of N'^ with H vanishes. Then the solution for the expectation values can be written as 

{L\t)) = {L')t=o + l{{L^)t=o + (iV')t=o)(exp(4|a|^M) - 1), (5.6) 

{N'{t)) = {N%^o + l{{L')t=o + {N')t=o)ieMMc^M - !)• (5-7) 

This result shows that the expectation values of L'^{t) and N'^{t) increase exponentially 
with time. If the Hamiltonian H is a possitive function of and AT^ alone, also the energy 
would grow in time. In this case the master equation describes a physical system coupled 
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to a reservoir of infinite temperature which increases the energy and angular momentum 
exponentially. However, such a system is not applicable for a phenomenological description 
of nuclear processes. But it is thinkable that realistic applications in solid-state physics 
can be found. 

Let us assume that the system has the smallest possible value of L^, namely {L'^)t=o = 
0. Then the increase of and N'^ is given according to Eqs. (5.6) and (5.7): 

{L\t)) = ^{N')t=o{exp{4\a\'m) - 1), {N\t)) = ^{N%^o{exp{4\a\^M) + 1). 

From these equations it is interesting to note that both {Lp'{t)) and {N'^{t)) are directly 
proportional to {N'^)t=o- Both expectation values are in general related by the equation 



{L'{t))-{L\t)) = {N'{t))-{N\t))t=o. 

(3) In the case Vj = apj + bqj, we obtain 

1 ^ 1 ^ 

L{A) = ^[H, A] + ^Dpp ^[g,-, [A, g,]] + ^D,, ^[pi, [APj]] 



h 



1 ^ iX ^ 

-^D,, [A,p,]] + b„ [A, g,]]) + ^([g,, [A,p,] + [A,p,]g, -p,[A, g,] - [A, g,]p,). 

j=i j=i 

(5.8) 

In this equation we used the following abbreviations: 

^qq = ^^*"' ^PP = ^^*^' ^pq = -^I^e(a*6), A = -Im(a*6). 

The quantities Dqq, Dpp and Dpq are denoted as the diffusion coefficients and A as the fric- 
tion constant. With the complex numbers a and b we can form the following determinant 
which is zero: 







h 


a*a 


a*b 




2 


b*a 


b*b 





Dqq -Dpq - iTlX/2 



-Dpq + ifiX/2 j^pp 

This determinant can be used to write down an equality serving as a constraint for the 
phenomenological parameters: 



DqqDpp - Dl^ = n^Xy4. 



(5.9) 



60 



If we set A = in Eq. (5.8), we find 

L(L2) = 4{Dppq^ + D,y - Dpgidp + pq)) - 4X{L'' + ^n^), (5.10) 

with q'^ = qf + + Qs^p'^ = Pi + pI + pI ^^'^ = QiPi + Q2P2 + QsPs- Since Eq. (5.10) 
contains other operators in addition to we have to calculate the corresponding equations 
of motion of these opertors, too. We obtain 

dq'^it) i 



dt n 



[H,q'{t)]-2Xq'{t) + 6D,,, (5.11) 



=-[H,p\t)]-2Xp'{t) + 6Dpp, (5.12) 

Let us first assume that the Hamiltonian depends on only, e.g., H = where 
is the moment of inertia. Then all commutators with H vanish in Eqs. (5.11)-(5.13). 
As function of time the expectation values of q^^p^ and qp + pq drop exponentially down 
(~ exp(— 2At)). For large times they approach the final values 

{q^)t^oo = 3-Dqq/A, {p^)t^oo = ^Dpp/X, ^(qp + pq)t^oo = ^Dpg/X. 

If we insert these values into the equation for d{L'^{t))/dt using Eq. (5.10), we find for 
large times 

(L2)t^oo = 0. 
Solving the expectation values we obtain 

{q''{t)) = Aexp(-2At) + 3DqJX, 

{p^{t)) = Sexp(-2At) + SDpp/X, 

^((qp + pq)t) = Cexp(-2At) + SDpJX, 

o 2 
{L^{t)) = L'exp(-4At) + ^{ADpp + BDqq - 2CDpq) exp(-2At). 

A 

The parameters A, B, C and D have to be chosen by initial conditions. A second example 
is the choice of the Hamiltonian of a spherical oscillator: 

" = 1/^ f 
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If we calculate the expectation value of L'^{t) for large times using the equations (5.10)- 
(5.13) we obtain 

<^'>'— = - rn^'D^^f+^'Di,}. (5.14) 

This relation can be rewritten by using the expectation value of the energy at large times: 

(iy),^oo = ^i^+mco^D,,). (5.15) 
Inserting (5.15) into (5.14) and applying (5.9) we obtain 

If the diffusion coefficients are chosen in accordance with Eq. (5.9), 

^pp ~ 2muj ' ~ ^pq — 

we reach the lowest value of {H)t^^ = (3/2)^a; and {L^)t-^oo = 0. 

In simple examples we have found operators V} needed to open a system in which 
angular momentum is damped. The difficulty in choosing these operators is the fact 
that in general new and usually complex operators appear on the right-hand sides of the 
Lindblad equations after evaluating the commutators. For these new operators one has to 
determine their master equations. Only if the number of coupled master equations remains 
small, one can hope to find analytical solutions for the damping processes. 

As shown in Eqs. (5.10)-(5.14) the damping depends on the parameter A which results 
from a product of two free parameters. If only a single free parameter is contained in the 
operators Vj, as for example given in the case (2), we obtain no damping of the square of 
angular momentum in general but an exponential increase. Systems with an energy and 
angular momentum increasing in time are not suitable for application in nuclear physiscs. 

A damping of angular momentum may be easier reached, if the rotational symmetry 
of the dissipative operators is given up. But the equations which have to be solved get 
more complex and cannot analytically be evaluated without approximations [88] . 

One can apply the damped three-dimensional harmonic oscillator presented here to 
heavy-ion collisions. For deep inelastic reactions of not too heavy nuclei one may approx- 
imate the internuclear potential near the barrier by a reverse parabola: 

711 

V{r) = VB--n\r-rBf. (5.16) 
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Here, r is the internuclear distance and tb the location of the barrier with height Vb- 
Equation (5.16) means a replacement of the frequency u> by ik. A damping of angular 
momentum of the relative motion can be gained if the friction constant A is larger than k. 
6. Quantum tunneling in open systems 

Tunneling is a remarkable, essentially quantum phenomenon, consisting in the prop- 
agation of a particle through a classically impenetrable potential barrier. Standing at the 
basis of important applications in electronics, chemistry and nuclear physics, it has con- 
tinuously been investigated from the beginning of the quantum mechanics [27], [89], [90], 
especially due to the difficulties raised by the effects of the dissipation [91-94], always 
present in practical cases. 

In this Section, the master equation (3.1) is used for the study of dissipative phenom- 
ena on tunneling processes. We consider a fission potential barrier U{q) [26], of the form 
represented in Fig. 4, where the first well corresponds to a compound nucleus and the sec- 
ond well to a fission channel. In order to generalize Gamow's formula of the tunneling rate, 
for the case when a dissipative environment in present we use the theory of perturbations. 
With the potential U{q), the Hamiltonian is 

H=^ + U(,) (6.1) 
and the master equation (3.1) becomes 

| = 4[J?.P]+AW, (6.2) 

where 

Hp) = -^i[Q,pp + pp] - [p,qp + pq]) 

-^Dgglp, Ip, p]] - -^Dpplq, [q, p]] + ^Dpq{[q, [p, p]] + [p, [g, p]]). (6.3) 

Using the commutation relation [g,p] = i^, from Eqs. (6.2)-(6.3) we obtain the 
following evolution equations for the expectation values of the canonical coordinates q and 
p: 



d<q> , <p> 
; + A < g >= — 
at m 

d<P> ^ dU ,^ ^, 

_^ + A<p>=<--> (6.4) 
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and for their variances: 



-^f^pp + 2A(jpp — 2(7p[7 = 2Dpp, (6.5) 

d 1 
do Tn 



where 



crqq^< q > - < q > , app ^< p > - <p> , aqq^< ^ > - < q><p> (6.6) 

and 

1 dC/ dU dU dU dU 

(^pu = -77 < -rP + p-r > + < -r ><P>^ (^qu = - < -i-q > + < — ><q> . 

2 dq dq dq dq dq 

(6.7) 

In this way we see that the quantum master equation (6.2)-(6.3) leads to the Newtonian 
equations (6.4) with additional friction terms and to the similar equations (6.5) for the 
variances, with Dgq, Dpp, Dpq as zero-point values and o'pq,o'pu,o'qU ^ generalized forces. 

Quantum tunneling is described as a transition from some "localized state" of the first 
well to some "localized state" of the second well. As we shall see later, this can be done by- 
defining a Hamiltonian Hq with approximate solutions in the first and in the second well, 
which disregard the tails corresponding to the other regions. The difference V between 
the two Hamiltonians defines a tunneling operator. Consequently, the Hamiltonian of the 
system takes the form: 

H = Ho + V, (6.8) 

where Hq describes the localized states and V the transitions between them. 

From Eqs. (6.4), with the expression U {q) — —Lolq'^/2, m = 1 for the potential barrier, 
one obtains: 

d<p> ^ -X<p> +0J^ <q> 
d< q> -X<q> + <p>' 

which can be compared with the similar equation standing at the basis of the Dekker's 

theory of the quantum tunneling with dissipation [28] , where the majority of the previous 

results are incorporated: 

^<P> =-2X+^^l^. (6.10) 



d<q> <p> 
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This equation has the solutions 

< p >i= i^ul + X^-X)<q>, < p >2= -{^(^l + A2 + A) < 9 > . (6.11) 

In this theory, the tunnehng rate is considered of the form: 

T{E) = Kv{E)P{E), (6.12) 

where P{E) is the penetrabihty and the coefficient k describing the effect of the dissipation 
is defined by the expression: 

«=^4^. (6.13) 

< p(0) > ^ ' 

With the first solution (6.11) describing the particle passing through the barrier, the ex- 
pression (6.13) takes the explicit form: 



« = ,/l+^-A, (6.14) 

which, in accordance to all previous results shows that due to the dissipation, the tunneling 
rate decreases. At the same time, with our expression (6.9) we obtain k = 1. We conclude 
that the result (6.14) is a consequence of the particular way of introducing the dissipation 
in the equation (6.10), which is only empirically justified. When one considers the more 
accurate equations (6.4), resulting from the quantum theory of open systems and which can 
be easily understood from a physical point of view, the effect described by the expression 
(6.14) disappears. As we shall see in the following, the dissipative environment generates 
new transitions increasing the barrier penetrability. 

By considering a basis of localized states, l^^o > for the ground state of the first well 
(the compound nucleus) and \^i>, « = 1, 2, ... for the second well (the reaction channel) 
the transition operator V takes the form: 



(6.15) 





( ° 


Vol 


V02. 






Vio 





0.. 





V = 







0.. 












0.. 


/ 



Vio = Voi, z = l,2, ...,n, 
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where the matrix elements have the approximate expressions [26]: 

Voi = = 2;^[*i(r-o)*i(ro) - *o(r-o)*;(ro)]. (6.16) 

We consider that the operators V and A(p) are smaU perturbations of the Hamiltonian 
Hq defining the basis of the locahzed states. With the notations: 

p(t) = e*"°'p{t)e-*^°\ V'{t) = e*^°'V{t)e-^^°\ A {p {ty,t) = e^^°*A(p)e-^^o* 

(6.17) 

in the interaction picture the quantum master equation (6.2) becomes: 

^ = -l[v'{t),p{t)]+A{p{ty,t). (6.18) 

Considering the zero-order density operator p = |0 >< 0| , for the initial state the 
density operator in the second order approximation with respect to the operators V and 
A, can be written in the form: 

p'(t) = p'(0)+pW(t)+p(2)(t). (6.19) 
Introducing this expression in Eq. (6.18) and integrating, we obtain: 

p{t) = p'(0) + py{t) + pj,{t) + p'yyit) + Pyj,{t) + p^yit) + P^J,{t), (6.20) 



where: 



Py{t) = ^y^dt[V'{t),p('^] 

p'jyit)^ fdtA{p^'^;t), 
Jo 



Jo 



Pvv{t) = -^ f dt f dt"[v'(t'),[v'(t"),p(o)]], 







P'VA = jy^' ) ' ^ 

Pav(^) = -^£^^'£ dtA{[v\t),p^'\t), 



t rt 

'(0). "- 



Paa(^)= dt dt A(A(p^^^t )-t). (6.21) 

'0 -'0 
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In principle we have to calculate the diagonal matrix elements Pii{t) in the representation 
of the localized states \ipi{t) >= exp{—iHQt/?i)\i >,i = 1,2,... In this case, Pii{t) =< 
il^i{t)\p\il^i{t) >=< i\p' {t)\i > and the transition rates have the expressiones Pii{t)/t . The 
functions pii{t)/t describe the dependence of the transition probabilities on the frequency 
Ui = {Ei — Eq)/% due to the energy uncertainty corresponding to the time t. In this 
way we can see the influence of the openess on the tunneling spectrum and on the total 
tunneling rate F = YliiPai^)/^- Using the expressions (6.3) and (6.15) of the openess and 
tunneling operators, neglecting the sums of the rapidly varying terms e*'^^* and making 
the integrals (6.21), from the expression (6.20) we obtain: 

Pii{t) = pfS) + Pu{t) + pUt) + P?i + Pii +Pti+ Pfiit) + Pfiit) + Pfiit) (6.22) 



with 



Piiit) = ni 



2 sin^(u;jt/2) 



iO' 



(6.23) 



pfiit) 



XCio ( XCoi + 2DggUoi - 2^«^oi) - 2^Qio{Coi + 2voi) 

Pii{t) — AOioiiOi JT, 1, 



n / \ c\/ ^vx> \ t . sin uj{t . 

Piiit) = -^{DqqUQQ + -^WQQ)CiQVLQi — { COSCUit), 

Tl LOi UJit 

pZ{t) = -2^XQo{Coi + 2^oi)-(l - ^), 
pi{t) = (^-2A?iosoi + 2Dggsl + 2^9'o) t, 

pfiit) = (^-XqiOSoi + Dqqsli + ^QiO^ Xt'^ , 



pfiit) ^Cor{2 



Pii — ^iQi 
f^Oi ^(C'oi + 2vQi) 



sin^(c^it/2) 



depending on the following overlap integrals or matrix elements: 



Coi = Cio = / *o(g)*i(g)c?g, ^^o* = ^i^ = 



Vio 

h ' 



(6.24) 
(6.25) 

(6.26) 

(6.27) 

(6.28) 

(6.29) 
(6.30) 



(6.31) 
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QzO = I '^Q{q)q'^^{q)dq, So^ = j ^'o(9)^^*^(9) 



uoi = - J -^-^^^^ ^00 = j (-^) 

V0i = j ^o{q)qd^i{q),Woi = J ^o{q)q^^i{q)dq,Woo = J ^o{q)q^^o{q)dq. (6.32) 

First of all, wc notice that although the matrix elements given by the expressions (6.22)- 
(6.32) depend on the very small arbitrary distance Sco between the energy levels in the 
fission channel, the sum: 

^Pii(t)~ J pii{uji;t)-^ (6.33) 

i 

of the probabilities on every span of energy do not depend on it, because every term 
contains two times the wave functions ^'i(g) which, according to the expression (6.24) is 
proportional to y/Suj . 

Secondly, one notice that, using the relation (6.33), the first term (6.23) leads to the 
Fermi's golden rule (6.32): 

T^ = ypM = 2n$^, (6.34) 

i 

where we have neglected the dependence on i of the matrix element Qio, in the very narrow 
span of energy Aa; = 1 /t. This term is symmetric with the energy Tuvi = Ei — Eq and has 
a very small width of the spectrum, representing the energy uncertainty for the time t. 

The second term (6.24), having the same time dependence as (6.23) is merely a cor- 
rection to the Gamow's term, due to the openess parameters. 

The third term (6.25), being asymmetric with uJi, describes a tunneling energy shift 
proportional with the coefficient A. Because the matrix element Qio is negative, this term 
increases the transition rates for energies Ei smaller than the ground energy of the nucleus 
Eo{(jiJi < 0), while for E^ > Eq the transition rates are decreased. 

The fourth term (6.26), rapidly varying with a;,, describes transition rate modifications 
due to the diffusion coefficients Dpp and Dqq which, on the other hand, do not change 
the energy expectation values, the new introduced positive terms being compensated by 
negative ones in the summation. 

The fifth term (6.27), depending on -Dpg, introduces an energy shift similar to that 
introduced by A. Because we have previously shown that at thermal equilibrium Dpq = 
[17], this term describes nonequilibrium processes. 
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The sixth term (6.28) leads to a transition rate pfi{t)/t also not depending on time, 
but with a very large spectrum, having a dependence on Ui much weaker than that of the 
functions multiplying the matrix elements in the expressions (6.23)-(6.27). 

The seventh term (6.29) leads to a transition rate pf^{t)/t also with the very large 
spectrum given by the matrix elements qio,Soi, but proportional with time. From the 
fundamental constraints (3.4) we find that both terms (6.28) and (6.29) are positive. 

The terms (6.30) and (6.31) lead to a time-independent desintegration probability due 
to the overlap of the initial state |0 > with the states \i >. 

For a long interval of time, the evolution of the system can be described by the 
expressions (6.22)-(6.31). In this case, the validity condition of the perturbation theory 
Piiit) « Poo(O) holds also for large values of time, when J2i Puif) » Poo(O), the number 
of energy levels i being very large. Making the summation of the expressions (6.22)-(6.31), 
by using the formula (6.33), one obtains an expression of the form 

Y,Piii^)-X + ^ot + Tit', (6.35) 

i 

where: 

Tot = E [pfiit) + pf,(t) + pUt) + Puit) + pZ it) + piit)] , 

i 

Tit' = E/^^^W'X = E + P^^^)] ■ (6.36) 

i i 

If N{t) is the number of nondesintegrated nuclei, we can consider the equality 

which leads to the following desintegration law : 

N{t) = iVoe-^-'^°*-^i*' = iV(0)e-^"*-^i*'. (6.38) 

In this expression appear two additional parameters: Fi describing the irreversibility of the 
desintegration process in a dissipative environment and x describing the nonorthogonality 
of the initial state |0 > of the nucleus with the fission channel states \i>. Estimating the 
time to = x/ro ~ mAq^/h ~ 10~^^s, where Aq is the width of the barrier, we find that x 
can be neglected. 
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In the tunneling spectrum described by the expressions (6.22)-(6.31) we distinguish 
the Gamow's spectral line: 

/ cjit) = , ^ ; ' 6.39) 
represented in Fig. 5 and additional lines due to the openess, described by other functions: 



g'^iujit) = — - cos UJit , (6.41) 



1 sin UJit 


— cos UJit 


UJit \ 


sin ujit\ 


UJit J 



(6.42) 



represented in Figs. 6, 7 and respectively 8. In this case, the transition rate spectrum can 
be written in the form: 

Tiuji) = P^ = Wf^tF{ujit) + ^{uji){l + ^t\ (6.43) 

where: 

Wl = nl + XCio (^XCoi + 2D,,uoi - 2^«^0i^ - 2^nio{Coi + 2voi) ~ Q^, (6.44) 

^{u^) = -2\q,oSo^ + 2Dqqsli + 2^q% > 0, (6.45) 

F{ujit) = fiujit) + rj^g\ujit) + v'^Q^'M + V^'g^'M, (6.46) 

the quantities 



^ = -X^^ -X-f > 0, 



D o / n I ^PP \ Cio^Oi o / n i ^PP \ \ n 

r] = -2 yDggUoo + ^2~'^00 J ^2 - ~^ yDqqUoo + -^^00 J ^ > 0, 

^ = ^D^^ CMCo^ + 2vo^) ^ Q 

being practically independent, positive parameters, describing the four processes of the 
interaction of the system with the environment: the friction, the diffusion of q, the diffusion 
of p and respectively the non-equilibrium. The function (f{uji) is also positive due to the 
fundamental constraints. 
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7. Open quantum systems and the atom-field interaction 

— * 

We consider a linearly polarized single mode electromagnetic field E, of the frequency 
u! and the wavevector ki propagating through an absorbing medium of two level atoms 
with the transition frequency ujq, the electric dipole ^ and the density A^; every atom n is 
described by the Pauli operators u^, (j^ , cr^. Due to the interaction between the radiation 
and the atoms, higher order harmonics Ei, of the frequency i/u and wavevector k^, are 
generated. For the electric dipole interaction, the Hamiltonian of the system is of the form 

H = + - ^{pI + v'^u'^ql) ^- ii^al{vuq^ sinker n^-Pu cos kj,r„), (7.1) 



where Ji, = /i/^/ioV", V is the volume of quantization, p,^, qi, are the canonical variables of 
the harmonic i/ , and fn the position vector of the atom n. In this case, we consider the 
operator Vj from the equation (3.1), of the form [29]: 

v n 

With the expressions (7.1) and (7.2) we obtain the master equation: 

^ = ~^ [^0, p] - ^[{^^Q<^ sin k^r^i + Pu cos k„rn)cr^, p] 



nv 



^ -ill n n 



{[q., [p., P]] + [p., [q., P]])} + [<' ^]] - KzK^ 



-ALK, [<, p]] + r:,([^:, K' P]] + K, p]]) + r-,([<, [a-, p]] + K, [a-, p]]) 

+rL(K, K, P]] + K, P]]) - liD^iWl a^P + pa^] - <P + P<]) 
+^^(k^ <P + P^^] - K, <P + P<]) + D-iK, a-p + pa^] - a-p + pa^])]}, (7.3) 

where Hq is the Hamiltonian of the system without interaction, and the phenomenological 
parameters are defined by the expression: 

3 3 
71 



j 3 
j j 3 

3 j 

^yz ~ '^^{Bjn^jn + ^jnCj^), — ^^(-^jn^jn ~ BjnC'jn): 

3 3 

^zx ~ y~l(^.7*n^jn + CjnA*^), Dy = — — ^(Cj„^j„ — CjnAj^). (7.4) 

i i 
Separating the field from the atomic observables by the "mean-field approximation", for 
the atomic expectation values we obtain the equations 

— < (7a; > +7^ <crx> -{(^o - s) <ay> +71 < cr^ > -Di = 0, 

^<cry> +{(jJo - s) <ax> +7^ < CTy > +(72 - x) < cTz > -D2 = 0, 

j^<(Tz> +li<(yx> +(72 + X) + 7|| < > -^3 = 0, (7.5) 

where 

7^ = 4(Aa;j/ + A^cc), 7I = ^{^xy + Ay^), 7|| = 4(Azx + Ay^), 
s = 4r,;j„ 71 = 4r,a., 72 = 4rj,„ i^i = 4^^, L>2 = 4L>y, L>3 = 4L>^ (7.6) 

and 

v = 



is a normalized field variable. From these equations one can obtain the classical Bloch- 
Feynman equations when s = 71 = 72 = 0, -Di = -D2 = 0, 7j^ = 7j^ [29]. 

Having in view the equation (3.3), in the master equation (7.3), we distinguish two ad- 
ditional parts describing the interaction with the environment: one part for the harmonic 
oscillators coresponding to the electromagnetic field and another part for the ensemble of 
atoms. This latter part describes three processes of the atom interaction with the environ- 
ment: 1) the friction, described by the parameters A or 7j^, 7j^, 7|| ; 2) the diffusion, leading 
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to zero-point values of the observables, described by the parameters Df, 3) the couphng 
between the observables described by the parameters F or s, 71,72- These parameters 
satisfy fundamental constraints, similar to the expressions (3.4) and (3.24): 

Dqqi/ ^ 0, Dppu > 0, Dppi/Dqqu -^^pgi/ ^ Z , 



/^2^2 

■Dqqi/^ppu(t) ~\~ Dppi,^qqjj{t) IDpqj^^pqjyit) ^ 



AS,>0, A-,>0, A-,>0, 

An An pn^ ^ ^7")"'^ A*^ A*^ "pn^ ^ ^ T")'^^ A'^ A*^ ■pn'^ ^ ^ T")"-^ 

xy zx yz — ' yz xy zx — 4 2^ ' •^^c y« a:2/ — 4 ' 

AS,A^,(t) + Kx^^yyit) - 2V-,^-M >d:< a^(t) >, 

A^.A-Jt) + A-^A:,(t) - 2r:,A:,(t) > i?" < a-(t) >, 

A^,A^,(t) + A-,A-,(t) - 2r^,A^,(t) > < a^{t) >, (7.8) 

where App, A^g, A^^^;, Aj,^, A^^^ are the variances of the observables and the quantities 
Apq, Aj;y, Ayz, A^^ ^1X6 exprcsslous of the form: 

Apg =< ^ >-<p><q>. (7.9) 

Following the model of the geometrical representation originally adapted by Feynman, 
Vernon and Hellwarth [95], we consider the atomic observables with the "amplitudes" 
u, V, w, as components of the Bloch vector in a "rotating frame" of the frequency uj : 

< CTx >= u cos u>t — V sin ujt, < ay >= —u sin ujt — v cos u>t, < ag >= —w. (7-10) 

For a system of atoms with the density N, using (7.10) one obtains the following expres- 
sions of the macroscopic polarization S = N < > and of the population N = N < >: 

S = ^(Se"^"^* + 8*6^"^*), N = -Nw, (7.11) 

where 

S = N(u-ii;). (7.12) 
With the expressions (7.10)- (7.12), the equations (7.5) become: 

^ + (7i + A)S + 7aS*e2^'^* ^ ^ _ ^gic^t ^ ■^*g2»a;t)^ ^ j^^iu^t ^ 
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^ + (AT -Ns) = l [{ix* - l*e'^' + ixe-^'-')S - {ix - le^^' + ix* e^'-')S''] , (7.13) 
where 

11. = ill + 7l)/2, la = (7I - 7l)/2, A = ujo-s-uj, 

7 = 71+^72, D = N{Di+iD2), iVa^NDa/Ty- (7.14) 

By neglecting the rapidly varying terms, one obtains the classical optical Bloch equations 
in the "rotating wave approximation" : 

-^+l\l{N-Ns) = -{Sx*-S*x). (7.15) 
In order to include the effects of these terms we consider solutions of the form: 

E = Eo + -(Eie"*'^* + E*e^'^*) + -(Eae-^*'^* + E^e^^'^*) + 

S = So + liSie-'''* + S*e^'"*) + ^(Sse-^*'"* + S^e^*'^*) + 

AT = No + -(Nie"''^* + N*e*^*) + -(Nae-^^'^* + N*e^''^^) + (7.16) 

where Eq, Ei, E2, Sq, Si, S2, Nq, Ni, N2, ... are slowly varying functions which are 
called "amplitudes". For the amplitudes of zeroth and first order we get the following 
equations: 

^ + (7^ + zA)Si = -7iNi, 

^ + 711 (No - ATg) + 71S0 = ^(Sixt - Stxi), 

^ + 7IS0 + 7i(iVo - iVi) = 0, (7.17) 

where xo = A*Eo/^, xi = A^Ei/^, -/Vi = '^D^hx . 

In these equations the new parameter 71 7^ leads to a coupling of the polarization 
variable Si with the electric field variable Xi, i-©- to an interaction of the atom with the 
electric field. In this case one obtains the polarization equation 

dS, ^ f 7l-^7;(72-2xo)/7i ^ ^ 2iXi^o 

dt \ 1 - H72 - 2xo)/7i / 1 - H72 - 2xo)/7i 
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which can be compared with the conventional Bloch polarization equation (7.15). We see 
that this equation describes new physical phenomena. For instance, in the steady state, 
with the notations 

71 7± " 7±7|| 2 

N^ = - 5=^, £=^, (7.19) 



1 - ^ 7± X. 



s 



7_l7|| 

one obtains the polarization amplitude 

" 7[ ■ (i + r5)2 + (5-C)^ + (i + r5)kP" ^ ^ 

For the slowly varying amplitudes xi+ ^ind xi- of the field variables, 

Xl = Xl+e^''+Xl-e-^'^ (7.21) 

the Maxwell equations take the form [96]: 

c^Xi+ ^ ^/f^g 

ci2; ch ' 

Xo = -^So, (7.22) 
Soft 

where ^ = LOji/^SQ is the coupling coefl&cients and Si+, Si_ are the Fourier transforms of 
the polarization Si as a function of the coordinate z: 

Si+ ^^j^ Si{z)e-^''^dz, Si_ Si{z)e"''dz. (7.23) 

Eqs. (7.22) with the expression (7.20) describe the propagation of an electromagnetic plane 
wave, taking into account the self-reflection effect [97] . Neglecting the counter- propagating 
wave due to the self-reflection, one obtains the absorption coefficient 

J_d\e[^^ l±Il (7 24) 

\s\ dz 2 (i + r5)2 + (5-c)2 + (i + r5)|£|2 ^'-^^^ 
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and the coresponding dephasing 

de s-c 



dz 1+T5 



a, (7.25) 



where ckq = 4:figN'^/ehjj_ . 

According to Eq. (7.24) the absorption coefficient a becomes negative when the 
electric field intensity |£p is sufficiently small for values of the atomic detuning 5 with 
1 + T5 < 0. In that case the electromagnetic wave is amplified, taking energy from the 
environment. At the same time, the diff'erence between the population of the ground state 
and the population of the excited state. 

No = z r-iT (7.26) 

^ + (<5-C)^+(i+r5)^ xi 

becomes larger than its equilibrium value N'^. When the ensemble of atoms has a smaller 
temperature than the environment, the energy passes from the environment to the atomic 
system and finally to the electromagnetic field. The absorption coefficient asymmetry with 
the atomic detuning S, can be observed in the experimental data of Sandle and Gallagher 
[27,98]. The laser radiation absorption by Na atoms in a buffer gas of Ar atoms has 
been studied by McCartan and Farr [99]. They varied the pressure of the buffer gas and 
measured the spectral line width w and the resonance frequency shift Aujq. Using the 
photoluminescence method [100], we obtain the following expression of the spectral line 
width 

= 27rAa;2-p^, (7.27) 
J aduj 

while the resonance frequency shift becomes 

Aa;o = Tie (7-28) 
Integrating the expression (7.24) one obtains 

= 271(1 + re). (7.29) 
By taking into account (7.28), this expression becomes 

W = 2AuJo ^ ^ 2Aa;o. (7.30) 
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Consequently the spectral line width W is proportional to the line frequency shift Aa;o 
which is in agreement with the experimental results of McCartan and Farr [99] . 

8. Summary 

The Lindblad theory provides a selfconsistent treatment of damping as a possible 
extension of quantum mechanics to open systems. In the present review first we studied 
the damped quantum oscillator by using the Schrodinger and Heisenberg representations. 
According to this theory we have calculated the damping of the expectation values of 
coordinate and momentum and the variances as functions of time. The resulting time 
dependence of the expectation values yields an exponential damping. Then we have shown 
that the quasiprobability distributions can be used to solve the problem of dissipation for 
the harmonic oscillator. From the master equation of the damped quantum oscillator we 
have derived the corresponding Fokker-Planck equations in the Glauber P, the antinormal 
ordering Q and the Wigner W representations and have made a comparative study of 
these quasiprobability distributions. The Fokker-Planck equations we obtained describe 
an Ornstein-Uhlenbeck process. We have proven that the variances found from the Fokker- 
Planck equations in these representations are the same. We have solved these equations 
in the steady state and showed that the Glauber P function (when it exists), the Q and 
the Wigner W functions are two-dimensional Gaussians with different widths. We have 
also calculated the time evolution of the density matrix. For this purpose we applied the 
method of the generating function of the density matrix. In this case the density matrix 
can be obtained by taking partial derivatives of the generating function. The generating 
function depends on a set of time-dependent coefficients which are calculated as solutions 
of linear differential equations of first order. Depending on the initial conditions for these 
coefficients, the density matrix evolves differently in time. When the asymptotic state 
is a Gibbs state in the case of a thermal bath, a Bose-Einstein distribution results as 
density matrix. Also for the case that the initial density matrix is chosen as a Glauber 
packet, a simple analytical expression for the density matrix has been derived. The density 
matrix can be used in various physical applications where a bosonic degree of freedom 
moving in a harmonic oscillator potential is damped. For example, one needs to determine 
nondiagonal transition elements of the density matrix, for an oscillator perturbed by a 
weak electromagnetic field in addition to its coupling to a heat bath. The density matrix 
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can also be derived from the solution of the Fokker-Planck equation for the coherent state 
representation. 

The theory was applied to the damping of the charge equilibration in deep inelastic 
collisions of heavy ions. The comparison of theoretical results with experimental data 
shows that the overdamped solution succeeded to describe these experimental data. 

We have also shown the Hamiltonian of the proton and neutron asymmetry degrees of 
freedom in deep inelastic collisions as an example for two coupled and damped oscillators. 
The usual limitation of the Lindblad theory is that the damping time is long compared 
with the characteristic time of the oscillators. This condition is not too well satisfied in 
deep inelastic collisions of nuclei, where the time scale is of the order of the relaxation time. 
Therefore, in these applications we consider the Lindblad theory as an axiomatic procedure 
for describing the dissipation processes and accept its parameters as free quantities, fitted 
to the experimental data. 

We obtained simple models for the dissipation of the angular momentum, by using 
the generators of the Lorentz group as openess operators of the system. For given values of 
the diffusion and friction coefficients we obtained an exponential damping of the angular 
momentum, in agreement with the situation encountered in heavy ion collisions. If the 
opening of the system fulfils rotational symmetry and the Hamiltonian is a positive function 
of and N'^, (square of generators of Lorentz group), then the Lindblad master equation 
describes a physical system coupled to a reservoir of infinite temperature, which increases 
the energy and angular momentum exponentially. 

We calculated the tunneling spectrum as a function of some barrier characteristics: the 
tunneling operator, the overlap integral, the transition elements of the coordinate and mo- 
mentum. Besides Gamow's tunneling process with energy conservation, additional terms 
with energy transfer to or form environment are obtained. For low values of temperature, 
the whole spectrum is situated at energies smaller than the Q-value; for higher enough 
values of temperature, transitions at energies higher than the Q-value appear. Generally, 
we shown that dissipation stimulates the tunneling process and leads to an exotic decay 
law. 

Using Lindblad's theory of open systems, for the resonant atom-field interaction, we 
found new optical equations. We showned that the interaction with the environment 
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consists not only in a decay, but also in a coupling of the atomic observables. This coupling 
has experimental evidence in the absobtion spectrum of the laser radiation and in optical 
bistability. 

At the same time, this phenomenon leads to a remarkable physical effect of energy 
transfer from the dissipative environment, which is cooling, to the coherent electromagnetic 
field propagating through it. 

Recently we assist to a revival of interest in quantum brownian motion as a paradigm of 
quantum open systems. There are many motivations. The possibility of preparing systems 
in macroscopic quantum states led to the problems of dissipation in tunneling and of loss 
of quantum coherence (decoherence) . These problems are intimately related to the issue of 
quantum-to-classical transition. All of them point the necessity of a better understanding 
of open quantum systems and all requires the extension of the model of quantum brownian 
motion. Our results allow such extensions and also explain some earlier observations. For 
a first comment of the result {3A),m) see Sec. 4.2 of Ref. [101]. The fact that the 
Schrodinger equation corresponding to the Lindblad equation is nonlinear was obtained in 
Ref. [17] and was quoted in Ref. [102]. The approach used in Ref. [17] to generate the 
Dekker master equation was applied in Ref. [93] to generate a family of master equations 
for local quantum dissipation. New developments of this kind were also obtained in Ref. 
[103]. 
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